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PASCO: STRUCTURAL PANEL ANALYSIS AND SIZING 

CODE , USER'S MANUAL 

Melvin S. Anderson, W. Jefferson Stroud, 

Barbara J. Durling, and Katherine W. Hennessy 

SUMMARY 

A computer code denoted PASCO is described. Buckling and 
vibration analyses are carried out with a linked-plate analysis 
computer code denoted VIPASA, which is included in PASCO. Sizing 
is based on nonlinear mathematical programming techniques and 
employs a computer code denoted CONMIN, also included in PASCO. 
Design requirements considered are initial buckling, material 
strength, stiffness, and vibration frequency. The report serves 
as a user's manual for PASCO. 

INTRODUCTION 

This report serves as a user's manual for a computer code 
denoted PASCO > Panel Analysis and Sizing COde. The code can be 
used to analyze and size prismatic structures having an arbitrary 
cross section. Its primary focus is stiffened panels made of 
laminated orthotropic materials. This tool is of particular value 
in analyzing and sizing filamentary composite structures. 

The capability of PASCO and the approach used in the 
structural analysis and sizing are described in reference 1. When 
used in the analysis mode, PASCO calculates laminate stiffnesses, 
lamina stresses and strains (including the effect of temperature 



and panel bending), buckling loads, vibration frequencies, and 
overall panel stiffness. When use' in the sizing mode, PASCO 
adjusts sizing variables to provide a low-mass panel design that 
carries a set of specified loadings without exceeding buckling or 
material strength allowables and that meets other design require- 
ments such as upper and lower bounds on sizing variables, upper 
and lower bounds on overall bending, extensional and shear stiff- 
nesses, an i lower bounds on vibration frequencies. 

The report begins with a discussion of structural modeling 
with PASCO. Plotting is then discussed briefly. Computer program 
definitions (including input and output) and organization are then 
presented. Finally, several illustrative examples are presented. 

MODELING WITH PASCO 

The modeling approach used in PASCO is basically the same as 
that used in VIPASA (refs. 2 to 4) . Although emphasis is placed 
on flat panels having several identical bays the only restriction 
on configuration is that the structure is assumed to be prismatic 
(fig. 1) . In addition, it is assumed that loads and temperatures 
do not vary along the length of the panel. Input is in NAMELIST 
format . 

The general procedure for setting up a problem is as follows: 

(1) Identify a repeating element of the structure. 

(2) Associate with each flat plate element a unique nvunber 

I and corresponding width B(I). (Curved segments may be modeled 
with several flat elements.) 
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Figure 1.- Typical prismatic structures that can be 
analyzed and sized with PASCO. 

(3) Associate with each unique ply a nvimber J. A ply is 

described by the following four input quantities: thickness, 

T(J); ply angle, THET(J); material, MAT(J,L); and temperature, 
TEM(J,L) . (Temperature and material properties may be different 
for each load case L.) 

(4) Identify each unique combination of plies that form a 
wall. The stacking sequence for a wall is prescribed with the 
input quantity KWALL(K,IW). 

(5) Specify wall for each plate element by using input 
quantity IWALL. 

(6) Specify properties of each material in namelist MATEk. 

(7) Create the geometry of the repeating element through 
HCARD input. 

(8) Specify number of bays with the input quantity NOBAY. 

(9) Define the loading on the panel: NX(L), NXY (L) , NY(L), 

NX(L) , PRESS (L) . 
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(10) Define the overall bow on the panel: ECC(L). 

(11) Define the frequency requirement: FR£Q(L). 

(12) Sizing variables are B(I), T(J), and THET(J) . Specify 
active sizing variables by inputting B, T, and THET as negative 
numbers . 

(13) Specify number of sizing cycles with the input quantity 
MAXJJJ . (MAXJJJ = 0 gives analysis only.) 

The above constitutes the minimum input required to handle a large 
number of cases. 

A complete set of input data and definitions is given in the 
section entitled COMPUTER PROGRAM DEFINITIONS AND ORGANIZATION. 

Much of the input is optional to give more flexibility to the user. 
Input that is essential is identified with an asterisk on the left 
side of the page in the section entitled Input. It is recommended 
that in the early stages of using the program, only essential 
input data be used to run problems. Such an approach requires 
only the thirteen basic steps listed above. 

Several aspects of modeling, including the input quantities 
KWALL, IWALL, HOARD, and ICARD, and two types of linking are 
discussed in the following sections. 

Basic Definitions 

Plate .- Each flat element defined by a width B(I) and 
associated wall IWALL (I) is defined as plate I. There are 
four degrees of freedom at each edge (three displacements u, v. 
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w, and a rotation 0 ) , so that an eight-by-eight stiffness matrix 
describes the structural characteristics of the plate. 

Substructures . - Substructuring is used to increase the 
efficiency of the analysis and to simplify the modeling of 
complicated configurations. Substructures are created with HCARD 
input, which is discussed in subsequent sections and illustrated 
with several examples in the final section of this report. Two 
types of substructures can be created: doubly-connected sub- 

structures and singly-connected substructures. 

• Double-connected substructure - A doubly-connected sub- 
structure is any assembly of plates and substructures with internal 
degrees of freedom removed and with only a beginning and final 
node. A plate is a special case of a doubly-connected substructure. 
Since a doubly-connected substructure has two nodes, its response 
is governed by an eight-by-eight stiffness matrix. Plates and 
doubly-connected substructures are numbered 1 to 120 and have an 
initial edge and a final edge that must be properly accounted for 
when assembling the final structure. The original orientation of 

a plate is horizontal witl. the initial edge on the left. 

• Singly-connected substructure - If boundary conditions are 
known at one edge of a substructure, those degrees of freedom may 
be removed and only a four-by-four stiffness matrix is required. 
Such a substructure is denoted a singly-connected substructure. 

It is essential to understand that the degrees of freedom on the 
initial edge are the ones removed and that only the final edge is 
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available for attachment to other nodes. Singly-connected 
substructures are identified by numbers 121 to 899. 

KWALL and IWALL Input 

The input quantity KWALij(K,IW) is a sequence of integers 
indicating the stacking sequence for wall number IW. The 
sequence begins with the outside layer. Since all walls in PASCO 
are symmetric, only one-half the wall is specified. Each integer 
is a ply number which defines a thickness, orientation angle, 
temperature, and material. A layer with a negative orientation 
angle is specified with a negative integer. For example: 

KWALL(1,3) = 5, -5, 2, 

means that wall number 3 has a stacking sequence of 

(THET(5), - THET(5), THET(2))^ 

s 

with associated thicknesses 
(T(5) , T(5) , T(2) 

D 

temperatures 

(TEM(5,L), TEM(5,L), TEM(2,L))^ 

O 

and materials 

(MAT(5,L), MAT(5,L), MAT (2,L))_ 
where L is the load case number. 
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The input quantity IWALL(I) is a sequence of integers indicat- 
ing the wall number for plate element I. For exeimple, assume 
that a stiffened panel is modeled with 7 plate elements having a 
total of 3 unique wall constructions. Assume that the pertinent 
input data are 

KWALL(1,1) » 1, -1, 2, A, 

KWALL(1,2) = 3, 6, 

KWALL(1,3) = 5, -5, 2, 

IWALL(I) = 1, 1, 2, 1, 3, 1, 2, 

The three KWALL vectors define the three unique walls. The IWALL 
vector means that plate elements 1, 2, 4, and 6 have wall 1; plate 
elements 3 and 7 have wall 2; and plate element 5 has wall 3. 

Specifying only the unique walls with KWALL and IWALL 
decreases storage requirements and run time. 

Boundary Conditions 

As explained in reference 1, various boundary conditions can 
be specified on the lateral edges of a stiffened panel, but 
boundary conditions cannot be specified on the ends of a panel. 

The simplest way to apply boundary conditions on the lateral 
edges of a stiffened panel is to use the input quantity IBC. The 
boundary conditions availedile with IBC are presented in table 1. 
When IBC is used, the same boundary conditions are applied to 
both the initial edge and the final edge of the panel. 
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TABLE 1.- BUCKLING BOUNDARY CONDITIONS ON LATERAL EDGES 
OF PANEL PRESCRIBED BY IBC 


IBC 

Boundary Condition 
Name 

Restraint on rotation ^ 

and displacement i 

0 

Free 

6, w, V, u are unrestrained 

1 

Simple Support 

6 and v are unrestrained, 
w = u = 0 (±909.0900) 

2 

Symmetry 

w and u are unrestrained, 
C = V = 0 (±990,9000) 

3 

Clamped 

6=w = v = u= 0 (±999.9900) 


If IBC input is inadequate, HCARO or ICARD input can be used 
to prescribe boundary conditions at any node. Boundary conditions 
are entered with two integers as ±9XX,XXOO where the four X's 
indicate restraint against (1) rotation 0, (2) displacement w, 

(3) displacement: v, and (4) displacement u, respectively. 

The displacements u, v, and v; are in the overall panel 
coordinate system shown in figure 2. An X = 0 means no restraint, 
and an X = 9 means an infinite restraint. See exeunples in 
Table 1. Without any specification the node is unrestrained, 
which corresponds to a free edge boundary condition. Xt is, 
therefore, not necessary to specify a free edge by making all X’s 
equal to 0. In fact, such a specification causes a program 
error. 

Examples presented in the sections entitled HCARD Input and 
ICARD Input illustrate how boundary conditions can be applied. 
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Z,w 

Figure 2.- Axes and displacements in overall panel 

coordinate system. 

HCARD Input 

Whereas KWALL and IWALL input are used to specify the wall 
construction for each plate element, HCARD input is used to 
manipulate these plate elements to create the geometry of the cross 
section. If ICARD input is used, all substructures used in 
ICARD input must be defined by HCARD input. 

The HCARD input is a continuous vector of integers that 
are separated into groups, or sequences, that correspond in meaning 
very closely to the original HCARD input to VIPASA (ref. 3) . In 
VIPASA, each sequence is on a separate card. 

With HCARD input, new plates or substructures can be defined 
from those previously generated by (1) rotating a previously 
defined plate or substructure, (2) offsetting nodes associated with 
initial and final edges of a previously defined nlate or substruc- 
ture, and (3) joining previously defined plates or substructures. 

Denote the numbers on an HCARD sequence by i^, -^ 2 * ~^3' 

. . . tin* Then 
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Number of additional numbers in the sequence 
Number of plate or substructure being generated, 
ij Number of previously defined plate or substructure 
that is being modified or joined to another plate 
or substructure. 

i^ n > 3. Additional instructions involving rotation 

of el<jnic;»5:s . offsets, joining a series of plate 
elements or substructures, and boundary conditions. 
These topics are discussed in the following 
sections . 

It is essential that the last HCARD sequence generate the geometry 
of the repeating element. 

Rotation .- If the orientation of a plate appearing in the 
final structure is not horizontal with initial edge on the left, 
it must be rotated to its final position. Rotation is accomplish- 
ed by the sequence 


HCARD — 4, 

where the initial integer 4 means there are 4 more integers 
making up the sequence, and i^ are all positive integers. A 
new plate or substructure number 12 is generated by rotating 
plate or substructure i^ through the angle which is defined by 
i^ and i^. 

If i^ IS greater than 50, 6 is given as in VIPASA by 

-^5 

6 » ±i^ X 10 ^ 


( 1 ) 


( 2 ) 
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For example, HOARD = 4, -5, 1, 90, 0, means rotate element 1 
by +90° and rename it element 5. This example is shown in 
figure 3. As another example, HOARD = 4, -6, 3, 450, -1, means 
rotate element 3 by +45° and rename it element 6 . 


Y 




Z 


HCARO • 4, -5, 1. 90, 0, 


Figure 3.- Rotation with HOARD input. Arrow points 
from initial edge of plate element to final edge. 

If i^ is not greater than 50, 0 is given by 


e * cos“^ ^ + 180°6 

±B(i5) 

-1 

where the angle defined by cos is between 0° and 180°, 

±B(i5) 

and 


6=0 f sign before i^ is positive 

6=1 if sign before i^ is negative. 

For exeuRple, HOARD = 4, -13, 10, -5, 8, means rotate element 10 
by the angle 



% 


9 = cos"^ “tiff 
and rename it element 13. 

Nodes . - The junction or attachment points of plate elements 
and substructures are called nodes. In addition, boundary condi- 
tions are applied at nodes, and plots of buckle mode shapes are 
based on deflections and rotations at nodes. When first defined, 
each nlate el nodes — one node at the initial edoe 

and one node at the final edqe of the plate element. These two 
nodes are indicated by the solid circular symbols at the edges of 
the plate elements in figure 3. Nodes also exist at the initial 
and final edges of a doubly-connected substructure and at the 
final edge of a singly-connected substructure. 

Nodes can be moved from the edges of plate elements and 
substructures, and they can be eliminated. Moving a node from 
the edges of plate elements and substructures is referred to as 
offsetting the nodes. Offsetting is discussed in the next 
section. Nodes are eliminated by making a plate or substructure 
part of a larger substructure. Eliminating nodes is discussed 
in a subsequent section entitled Joining plates or substructures. 

Offsets. - Conventional modeling, where plate element center- 
lines are connected, results in overlapping areas or missing 
areas. Offsetting nodes from the edges of plates provides a more 
accurate model of the intersection of two or more plates. Also, 
offsets allow unsymmetric laminates to be created by stacking 
symmetric laminates. See, for example, figure 7. 
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The nodes of a plate or substructure can be displaced 

from its end points to create a new plate I 2 by the following V 

sequence . 

HOARD = 6, -ij, ±i^, ±i^, ±ig, +i^, (4) 

The offsets in the y and z direction at the in;.tial edge are 
±T(ig), respectively, and for the final edge are 
tT(i6) ±T(i^), respectively. The sign convention is such 

that if the y-z coordinate system is placed at the center of the 
plate edge, the coordinates of the node are the proper values of 
the offsets. 

For example, 

HOARD = 6, -8, -1, -5, 0, -6, 7, 

creates a new plate 8 from old plate 1 with the offsets shown in 
figure 4. Note that the integer zero produces a zero offset. 

I ►Y »- —7 ■» • — -J— L 

U J ™ 

•T(6)-*l I*- ^ 

Z 

HCARD • 6. -8. -1. *5, 0. -6, 7, 

Figure 4.- Offsetting nodes with HCARD input. Arrow points from 
initial ecge of plate element to final edge. 
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The procedure used to offset the nodes at the edges of a 
doubly-connected substructure composed of many plate elements is 
the same as that used for the single plate element shown in 
figure 4. For a singly-connected plate or substructure, the only 
node is located at the final edge. Offsets for that node are 
defined by ±i^ and - 15 / not ±ig and -iy' although entries 
for ig and i^ must be included in the sequence. 

Joining plates or substructures .- Joining plates or sub- 
structures means taking previously defined plates and substructures 
and joining them to create the desired cross section. As explained 
earlier, plate elements and substructures are joined at the nodes 
associated with the edaes of the plate elements and substructures. 
For example, if plate elements and substructures are defined as 
shown in figure 5(a), they can be joined to foirm the cross section 
shown in Figure 5(b). Nodes are shown at the edges of each plate 
element and substructure. Arrows point from the initial edge to 
the final edae of each plate olement or substructure. 

Whereas HOARD input specifying rotations or offsets requires 
a specific format (eqs. (1) and (4)), HOARD input which specifies 
the joining of plate elements or substructures is relatively 
flexible in its format. The general format is 

HOARD ~i^,i. 2 ,i^,~^^,~i^,*** (5) 

where i 2 and i^ are positive. The rules for using HOARD input 
to join plates and substructures are stated briefly in the 
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(a) Component plate elements and substructures 


1 2 3 7 8 9 



(bl Components joined to form desired cross section 


Figure 5.- Plates and substructures joined to form 
cross section of prismatic structure, 

following paragraphs and are explained or illustrated in 

subsequent examples. 

If equation (5) is used to create a doubly-connected sub- 
structure, the integer ±2 is 1 to 120. If equation (5) is 
used to create a singly-connected substructure, the integer i2 
is 121 to 899. 

Only substructures having nodes that form a simple chain can 
be constructed with HOARD input. A simple chain is a simple path 
that connects all the nodes. 

All substructures have nodes temporarily ntjtinbered 1, 2, 

3, ..., N^. In the HOARD sequence given by equation (5), the 
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entry is the number of the first plate or doubly connected 

substructure that connects nodes 1 and 2. If additional plates or 
substructures attach at node 1 or if it is desired to apply a 
boundary condition at node 1, then these additional plates or sub- 
structures or these boundary conditions are all entered using 
negative signs — in any order. If the additional substructure is a 
doubly-connected substructure, it is mandatory that its geometry be 
such that its final node coincide with node 2 as defined by the 
entry i^. If the geometry is not compatible, the program still 
executes, but the results are incorrect. In this case, if the 
plotting option is exercised, a nonfatal diagnostic is printed. 
Boundary conditions are entered as -9XX,XX00 where X's indicate 
restraint against rotation 0 and displacements w, v, u, respectively. 
An X = 0 means no restraint and X = 9 is an infinite 
restraint. 

After all elements connecting nodes 1 and 2 and boundary 
conditions and singly-connected substructures at node 1 are 
accounted for, the next number in the sequence is positive, 
indicating a plate or substructure connecting nodes 2 and 3 . 

Again, additional negative numbers indicate other connections or 
boundary conditions at node 2. 

This process is repeated until the substructure is described. 
There is a limit of 20 integers allowed in any given sequence. 

If more integers are required, an intermediate substructure can be 
defined and used to start a new sequence. 
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If a singly-connected substructure or a boundary condition 
is at the last node (Nv) enter a zero then the substructure or 
boundary condition preceded by a negative sign. 

For singly-connected substructures, must be greater 

than one, and for doubly-connected substructures, must be 

greater than two. A diagnostic in the progrcun is printed if this 
rule is violated. 

Examples are now used to illustrate and explain the rules 
just discussed. In the examples presented in figure 6, four 
substructures are created from plate elements 1 to 6. Plate 
elements (or substructures) 5 and 6 were created earlier from 
plate element 1 by HOARD rotation similar to that illustrated in 
figure 3. The arrows on the plates and substructures point from 
the initial edge to the final edge. The nodes shown in the four 
examples in figure 6 indicate the nodes remaining after each 
substructure is created. Interior nodes are eliminated in HOARD 
modeling. 

In the first example in figure 6, HOARD = 2, 121, 5, 
plate element 5 is changed from a doubly-connected substructure 
to a singly-connected substructure d' loted 121. The bottom node is 
eliminated because the unattached f e of a singly-connected 
substructure must be the initial edge of that substructure. 

In the second example, HOARD « 4, 7, 6, 4, 5, the initial 
edge of plate element 4 is joined to the final edge of plate 
element 6, and the initial edge of plate element 5 is joined to the 
final edge of plate element 4. The resulting substructure 
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EXAMPLE 1 


EXAMPLE 2 


EXAMPLE 3 


HCARD ■ 2, 121, 5. 


HCARO ■ 4, 7, 6, 4. 5, 



HCARD - 4, 8, 1. 2. -121, 


121 


SINGLY-CONNECTED SUBSTRUCTURE 121 


DOUBLY-CONNECTED SUBSTRUCTURE 7 


DOUBLY-CONNECTED SUBSTRUCTURE 8 


EXAMPLE 4 




• 

HCARD • 5, 9, 1. 2, -7, 3, 

7 



DOUP'Y-CONNECTED SUBSTRUCTURE 4 


Figure 6.- Four examples of HCARD input used to join 
plate elements and substructures. 

is denoted 7. As shown in the figure, interior nodes are 

eliminated in substructure 7. Substructure 7 is a doubly- 

connected substructure. 

The third excimple, HCARD = 4, 8, 1, 2, -121, illustrates how 
two substructures can attach at the same node when one of the 
substructures is singly connected. The initial edoe of plate 
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element 2 is attached at the final edge of plate element 1 . To 
attach a second substructure, here denoted 121, at the final edge 
of plate element 1, enter the substructure number with a minus 
sign. To attach additional singly- or doubly-connected sub- 
structures at the final edge of plate element 1, enter each 
additional substructure number with a minus sign. 

The fourth example, HCARD =5, 9, 1, 2, -7, 3, illustrates 
how two doubly-connected substructures can attach at the same two 
nodes. The initial edge of plate element 2 (a doubly-connected 
substructure) is attached at the final edge of plate element 1. 

To attach a second doubly-connected substructure, here denoted 7, 
at the final edge of plate element 1, enter that substructure 
number with a minus sign. The minus sign before the 7 indicates 
that a second doubly-connected substructure, here denoted 7, is 
also attached at the final edge of plate element 1. In this 
example, the rule requiring a simple chain of nodes means that 
doubly-connected substructure 7 must reattach at the final edge of 
substructure 2. This is the reason that intermediate substructure 
7 was created from plate elements 4, 5, and 6. 

If a singly-connected substructure is at the last node, enter 
a zero then the substructure number preceded by a negative sign. 

A more complex example which illustrates this modeling rule is 
presented in figure 7. The nodes for plate element 22 have been 
offset by one-half the sum of the thicknesses of plate elements 13 
and 22. Although plate elements 13 and 22 must have symmetric 
leuninates, the combination of plate elements 13 and 22 need not 
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(a) PREVIOUSLY DEFINED PIATES AND SUBSTRUCTURES 




123 / 


2 

5 / 6 \ 13 



\ / 22 

121 


122 


(b) SUBSTRUCTURE CREATED WITH 

HCARD • B, 23. 2, -121, 5, -122, -123 , 6. -18, -21. 13, -22, 0, -124, 

Figure 7.- Example of HCARD input used to create 
complex substructure. 

produce a symmetric laminate. The two nodes shown in figure 7(b) 
indicate the nodes remain!. ig after doubly-connected substructure 23 
is created. 

Three examples illustrate how boundary conditions can be 
applied with HCARD input. In figure 8 the HCARD input 

HCARD * 8, 3, 1, -909, 0900, 2, 0, -909, 0900, 

places simple support boundary conditions at the node associated 
with the initial edge of plate element 1 and at the node associated 
with the final edge of plate element 2. The entry 0 after 
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hoard • 8. 3, I, -909, 0900, ?, 0, -909, 0900, 


SS .5 2, SS 

Figure 8.- Example of simple support boundary conditions 
applied with HOARD input. 

the 2 moves the second boundary condition to the final edoe 
of plate element 2. 

In figure 9, the HOARD input 

HOARD = 4, 121, 5 -900, 9900, 

4. 6, 1, 2, -121, 

places simple support boundary conditions at the tip of the blade. 
Note the difference between the -909, 0900 in figure 8 and the 
-900, 9900 in figure 9. The displacements w and v are in the 
z and y directions associated with the overall structure. To 
get a simple support boundary condition at the tip of the blade in 
figure 9, the v and w restraints had to be reversed from 
those specified in figure 8. 

In figure 10, the HOARD input 

HOARD = 2, 121, 5, 

7, 6, 1, 2, -999, 9900, 3, -121, 




o 

il 


HCARO • 4. 121, 5, -900, 9900, 
4, 6. 1, 2, -121, 



Figure 9.- Example of simple-support boundary condition 
applied at tip of blade with HOARD input. 



HCARO • 2 . 121, 5, 

7 . 6. 1, 2, -999, 9900, 3, -121, 



Figure 10.- Example of HOARD input used to apply 
boundary condition at interior node. 

places clamped boundary conditions at the interior node associated 

with the initial edge of plate element 2. All other edoes are 

free. 

Finally, there is a rule which requires that there be at 
least three nodes in a doubly-connected substructure and at 
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least two nodes in a singly-connected substructure prior to the 
elimination of any node for that HCARD sequence. For example, 

HOARD = 2, 3, 1, 

is not legitimate because there are only two nodes, and 

HCARD = 2, 122, 121, 

is not legitimate because singly-connected substructure 121 has 
only one node when it is used in this HCARD sequence. Also, 

HCARD = 3, 9, 2, -7, 

is not legitimate because substructures 2 and 7 connect the same 
two nodes. A diagnostic detects this input error. 

Reuse of plate and s ubstructure numbers .- In preparing HCARD 
input, the safest practice is to define new plates and sub- 
structures with numbers not previously used. However, numbers 
can be reused to save storage. Also, for very large problems, 
it might be necessary to reuse numbers to keep the element 
numbers within the maximum number of 120. Element numbers can be 
reused if certain rules are followed. When a new plate or sub- 
structure is defined with a number that has been previously used, 
the stiffness matrix for the first pl^te or substructure is 
destroyed so it is not available for further use. It is not 
permissible to define a rotated plate or substructure with the 
same number as the original plate. However, for offsets and 
joining together of plates, the new plate may have the same number 
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as one of the plates from which it is constr^icted. Finally, 
additional caution is required when using numbers that are not 
greater than the number of plate widths B(I) input. Original 
plate numbers may be redefined internally to increase the 
efficiency of the VIPASA analysis. If any of these numbers aie 
used for new substructures, a careful "'k of the HOARD data 
that VIPASA uses (which may be different from that input) should 
be made to see if the geometry is correctly defined. Plots 
showing the repeatinc element and total structure greatly 
facilitate this check. 


ICARD Input 

In most cases HOARD input is sufficient t" model structures 
being analyzed and/or sized with PASOO. lOARD input is used only 
when structures cannot be modeled with HOARD input or when a more 
detailed mode shape plot is desired. (Mode shape plots are based 
on displacement and rotation at nodes, and HOARD modeling 
eliminates interior nodes in a repeating element.) Even when 
lOARD input is not used, ICARD data are generated internally and 
are used by VIPASA. 

ICARD input is used to assemble the final structure from 
plates and substructures created by HOARD input. ICARD input is 
similar to HOARD input in that it consists of one continuous 
vector of integers grouped in sequences. Each sequence of 
integers is treated separate, y. 
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( 6 ) 


If a sequence of integers is denoted 


ICARD ij^, i2» i^f . . *^^+2' 

then 

ij^ Nvimber of additional numbers in the sequence. 

i 2 Node number. One sequence must be prepared for 

each node connected to a higher numbered node or 
having a prescribed boundary condition. 

i^ Node number (must be greater than i 2 ) that is 

connected to i 2 « 

Plate or substructure that connects the two nodes. 
Recall that each plate or substructure has an 
initial and final edge and it is essential that the 
initial edge be connected to a lower node number 
than the final edge. 

ij-i^ Additional pairs of numbers, the first being the 
number of the node which is connected to i^, and 
the second being the number of the plate or 
substructure making the connection. Note, if 
more than one element connects the two nodes, the 
higher numbered node number will be repeated. 

i^^^ Boundary conditions or the number of a singly- 

connected substructure at node i2 preceded by a 
minus sign. Notation for boundary conditions is 
the same as in the HCARD (±9XX,XX00) . 

i ^^2 Additional singly-connected substructures or a 

boundary condition, all entered as positive. 

An example which illustrates ICARD input is presented in 
figure 11. The ICARD input is given by 


ICARD = 5, 1, 2, 
3, 2, 3, 
3, 3, 4, 
3, 4, 5, 


1, -909, 0900, 

2 , 

3, 

4, 


3, 5, -909, 0900, 
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® 1 @ 2 ® 3' ® 4 ® 

ss* • ?SS 

ICARD = 5. 1, 2, 1, -909, (WOO, 

3, 2, 3, 2. 

3, 3, 4, 3, 

3, 4, 5, 4, 

3, 5, -909, 0900, 

Figure 11.- Example of ICARD input, including 
input prescribing boundary conditions. 

In the first sequence of integers, node 1 is connected to node 2 

with plate element 1, and a boundary condition is imposed at 

node 1. In the second sequence of integers, node 2 is connected 

to node 3 by plate element 2, In the fifth sequence of integers, 

a boundary condition is imposed at node 5. In each case, the 

plate elements could have been substructures created with HOARD 

input. 

For the analysis, the order of the ICARD input is immaterial; 
however, to obtain geometry and plotting, it is necessary that at 
least one node number used in any sequence after the first must 
have been used in an earlier sequence. 

When ICARD input is used, it is still necessary to define a 
repeating element with HOARD input in order to obtain the proper 
load distribution and to obtain necessary data for panel stiffness- 
es and plate centroid locations that are used in analyses involving 
bending loads. It is also necessary to define the boundary 
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conditions with ICARD input (±9XX,XX00) , because the IBC parameter 
is not used when ICARD input is used. 

The rules for ICARD input in PASCO are the same as the rules 
for ICARD input in VIPASA (refs. 2 to 4) with the additional 
capaibility of having to input only one repeating element of a 
repetitive structure. The input quantity ICREP takes the 
structure defined by ICARD input and repeats it ICREP times. 

(See section entitled NOSUB, NOBAY, and ICREP Input for specific 
input requirements.) 

An example which illustrates ICARD input combined with HCARD 
input, boundary conditions, and ICREP is presented in figure 12. 
The plate element components and the nodes are identified at the 
top of the figure. (Plate elements 2, 4, 7, 9, and 11 are 
assumed to have been rotated with previoiis HCARD input to 
produce plate elements 14, 15, 16, 17, and 18, respectively.) 

The repeating element, substructure 21, is first modeled with 
HCARD input. ICARD input is then used to provide detailed 
plotting. To retain nodes 3, 4, 6, 9, and 10, basic plate 
elements are used rather than substructures 19, 20 and 121 
created by HCARD input. Bounaary conditions are imposed with 
ICARD input at nodes 1 and 12 of the repeating element. The 
input ICREP * 2 creates the final structure shown at the bottom 
of the figure. In the final structure, there are two repeating 
elements, and the boundary conditions at node 12 are shifted 
to the final edge of the panel. In this example, all nodes are 
retained in the final structure. 



HCARD • 4, 19. M. 3. 15, 

4. 20. 17. 10. 18, 

2. 121. 16. 

10. 21, I, 5, -19, 6. 8, -121, 12, -20, 13, 

ICARD • 5. I. 2. I. -909, 0900, 

5. 2, 5, 5, 3, 14, 

3. 3, 4, 3, 

3, 4, 5. 15. 

3. 5, 7, 6. 

3. 6. 7. 16. 

3, 7, 8, 8, 

5, 8. II. 12, 9. 17, 

3, C , 10. 

3, 10, li, 18, 

3. II, 12. 13. 

3. 12, -909, 0900, 

ICREP ■ 2, 



FINAL STRUCTURE 


Figure 12.- Example of ICARD input, including HCARD and 
ICREP input and boundary conditions. 

Whereas only structures and substructures having nodes that 

form a simple chain can be modeled with HCARD input, that 

restriction is not imposed on structures modeled with ICARD input. 
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In the example presented in figure 12 , the nodal path 2 , 3, 4, 5 
and the nodal path 2 , 5 produce a double chain that cannot be 
modeled with HOARD input. To model the structure with HOARD input, 
intermediate substructure 19 was created. In doing so, nodes 3 
and 4 were eliminated. Nodes 3, 4, 6,- 9, and 10 can be retained 
only with lOARD input. Examples which illustrate complex modeling 
with lOARD input are presented in a subsequent section entitled 
EXAMPLES . 


NOSUB, NOBAY, and lOREP Input 

The input quantities NOSUB, NOBAY, and lOREP are used to 
replicate a repeating element generated by HOARD and/or lOARD 
input . 

HOARD requirements . - After creating a repeating element with 
HOARD input, an entire stiffened panel can be generated with the 
input NOSUB and NOBAY. The number of repeating elements in a 
major substructure is equal to , and the number of 

major substructures in the final structure is equal to NOBAY. 

The parameter lOREP is not used when the repeating element is 
created only with HOARD input. 

Assume that a repeating element with 1 stiffener has been 
defined with HOARD input. Also assume that a 16-stiffener panel 
is desired. The following combinations of NOSUB and NOBAY 
produce the 16-stiffener panel. 



NOSUB 


NOBAY 


1 16 

2 8 

3 4 

4 2 

5 1 

As a rule of thumb, a combination of NOSUB and NOBAY in which 
NOSUB and NOBAY are nearly equal is recommended for initial 
studies of a structure. The most efficient combination depends 
upon the structure. If a 15-stiffener panel is desired, the 
only combination is NOSUP = 1, NOBAY = 15. 

If NOSUB and NOBAY input are used and if boundary conditions 
are specified with HOARD ’nput (-9XX,XX00) , these boundary 
conditions are repeated on the corresponding node of each 
repeating element in the structure. An example which illustrates 
NOSUB and NOBAY input is shown in figure 13. 

ICARD requirements .- After creating a repeating element with 
ICARD input (which also requires HOARD input) , an entire stiffened 
panel caa be generated with the input quantities lOREP and NOBAY. 
The number of repeating elements in the final structure is equal 
to ICREP, which must equal NOBAY. The parameter NOSUB is not 
used with lOARD input. 

When using lOREP with boundary conditions, it is important 
to note that boundary conditions may be applied with two types 
of lOARD sequences, and that PASOO has been coded so that the type 
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HCARO ■ 4, 121, 5, -900, 9900, 
4, 6, 1, 2, -121, 


A 


NOSUB • 1 , 
NOBAY ■ 2, 
IBC- 1, 




■JS SS 


ki SS 


K SS 


FINAL STRUCTURE 


Figure 13.- Example showing how boundary conditions applied 
with HCARD input to repeating element are repeated at 
corresponding nodes of final structure and how boundary 
conditions applied with IBC input are applied only to the 
lateral edges of the panel. 

of ICARD sequence chosen will determine the location of the 
boundary condition. The two types of ICARD sequences are; 

(1) An ICARD sequence that contains node connections as 
well as boundary conditions. (For example, 5, 1, 2 , 
-909, 0900, in fig. 12) 

(2) An ICARD sequence that contains only boundary condi- 
tions - no node connections. (For example, 

3, 12, -909, 0900, in fig. 12) 

Boundary conditions prescribed with the first type of ICARD 
sequence are applied only to the first repeating element of a 
replicated structure. Boundary conditions prescribed with the 



second type of ICARD sequence are applied only to the last 
repeating element of the replicated structure. No boundary 
conditions are applied to the interior repeating elements. In 
the truss-core example presented in the section entitled EXAMPLES, 
boundary conditions are applied with ICARD input at more than one 
node at each edge of the panel. 

Linking 

Linking equations .- It is useful to be able to specify 
relationships involving all the B, T, and THET variables in 
order to define offsets, stiffener spacing, angles, and other 
quantities associated with panel geometry. A set of general 
linear equations relating all the variables can be input with the 
parameters AB(I,N), AT(J,N), ATHET(J,N), and AC(N). The linear 
equations in terms of these parameters are given by 

AB(I,1)*B(I) + AT(J,1) *T(J) + ATHET(J,1) *THET(J) *AC(1) 

AB(I,2)-B(I) + AT(J,2) *T(J) + ATHET (J, 2) -THET (J) =AC(2) 


AB(I,N)*B(I) + AT(J,N)*T(J) + ATHET ( J,N) *THET ( J) = AC(N) 

(7) 

Certain B, T, and THET variables in these linking 
equations are taken to be the dependent variables . The other 
B, T, and THET variables are taken to be the independent 
variables. Each dependent variable is input as 1.E30. A 
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consistent input then has the number of linking equations equal 
to the number of 1.E30 input variables. If the input is not 
consistent, a diagnostic message is printed. 

As an example, assume that the pertinent input is 

B(l) = -1.5, -1.5, 1.E30, 

AB(1,1) = 1., 2., 1., 

AC(1) = 6., 

The associated linking equation is 

B(l) + 2B{2) + B(3) =6 (8 

The negative signs on B(l) and B(2) indicate that they are 
sizing variables. Such negative signs are changed to positive 
signs for the linking equations. The initial value of B(3) is 
calculated by equation (8) to be 1.5. Because values of B(l) 
and B(2) may vary during sizing, subsequent values of B(3) 
may not be 1.5. (See discussion in example 3.) 

Automatic linking .- PASCO also contains an automatic linking 
feature that is applied only to plate element widths B(I) and 
is available when the parameter LINK is input as 1. Automatic 
linking is more convenient to use and saves storage and computer 
time relative to equivalent linking carried out with the linking 
equations discussed above. In general, automatic linking 
works in the followina way: After the linking equations 

(eq. (7)) are applied to the variables B(I), T(J), and THET(J) , 



plate elements B(I) are organized into linking groups having 
the same wall construction (same IWALL) and same width B. 

These linking groups are retained throughout the sizing. During 
sizing, the automatic linking feature equates the widths of all 
plate elements within a group to the width of the plate element 
having the lowest number in the group. 

For example, assume that the pertinent input data are given 
by 

LINK = 1, 

B(l) =1.5, 3. , 1.5, 1.5, 1.5, 

IWALL (1) = 1, 2, 1, 3, 1, 

Plate eleir-^nts 1, 3, and 5 form a linking group in which B(l) 
is the controlling width and B(3) and B(5) are equated to 
B(l) during sizing. Element width B(4) is not part of the 
linking group because IWALL(4) = 3 rather than 1. If a plate 
width is intended to be a controlling width, care should be taken 
to insure that it has the lowest plate nuniber in the group. For 
example, if the pertinent input were given by 

LINK = 1, 

B(l) =1.5, 3. , -1.5, 1.5, 1.5, 

IWALL (1) = 1, 2, 1, 3, 1, 

B(3) would always equal B(l) and would not act as a sizing 
variable. Caution — the default value for LINK is 1. 
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If a linking group contains more than one B that is a 
sizing variable, the sizing variable v;ith the lowest plate number 
is retained in the linking group, and all other sizing variables 
are removed from the linking group. For example, if the 
pertinent input were 

LINK = 1, 

B(l) = -1.5, 3. , -1.5, 1.5, 1.5, 

IWALL(l) = 1, 2, 1, 3, 1, 

the final linking group would include only B(l) and B(5), 
with B(l) as the controlling width. Both B(l) and B{3) 
would function as sizing variables. 

Linking equations (eq. (7)) are applied to the B(I) , T(J), 
and THET(J) before linking groups are identified. If linking 
equations cause a B(I) to be identified as part of a linking 
group, and if that B(I) is not the controlling width in that 
linking group, then automatic linking will override the linking 
equations in subsequent calculations of that B(I). For .example, 
assume the pertinent input were 

LINK = 1, 

B(l) = -1.5, 3., 1.E30, -1.5, 1.5, 

IWALL » 1, 2, 1, 3, 1, 

AB(1,1) » 0. , 0. , -1. , 1. , 
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The linking equation would be 


-B(3) +B(4) = 0 

With this input the user probably wanted B(3) to equal the 
sizing variable B(4) throughout the sizing. However, because 
the initial value of B(3) is calculated as 1.5, B(3) is 

identified as being part of the linking group composed of B(l), 
B(3) and B(5) with the sizing variable B(l) as the controlling 
width. If the user wishes to work the sizing problem with 
B(3) = B(4), the input for B(4) could, for example, be changed 
from -1.5 to -1.505. There are, of course, many other approaches. 

Load Distribution 

The longitudinal (N ) , transverse (N ) , and shear (N ) 

i ^i *^i 

loads in each plate element are automatically calculated based 

on the applied loads, the geometry and stiffness distribution of 

the cross section, and the modeling rules. The loads N and 

^i 

N in any plate element can also be prescribed using the input 
XYi 

quantities FNY and FNXY. 

Longitudinal load, N .- The load distribution N 
^ *i 

calculated by PASCO is the sum of a load distribution based on 

uniform longitudinal strain of all elements in the panel cross 

section and a load distribution caused by bending loads. The 

resulting load distribution cannct be overridden with input. 

Details of the stress analysis are presented in reference 1. 



Transverse load, N .-In general, the modeling rules are 
^ 

designed so that the full applied transverse load is carried 

by the skin of a stiffened panel, and no part of is carried by 

the stiffener elements. This objective is accomplished by the 

following two rules that are applied automatically in the program. 

(1) In a singly-connected substructure, such as a blade on a blade- 

stiffened panel, N is automatically set to zero. (2) If two 

^i 

or more doubly-connected substructures join the same two nodes, 

N of the first substructure mentioned on the appropriate HOARD 
^i 

sequence is automatically set equal to the full N , and N of 

y 

the second and subsequent substructures are automatically set equal 
to zero. For example, if HOARD input is 




HOARD = 5, 5, 1, 2, -3, 4, 


then substructures 1, 2, and 4 carry the full load, and 

substructure 3 carries no N load. The same structure could 

y 

also be modeled with the HOARD input 


HOARD = 5, 5, 1, 3, -2, 4, 


In this case, substructuie 3 carries the full load and 

substructure 2 carries no load. 

If these modeling rules do not provide the desired load 


distribution, the transverse loads in each plate element can be 
specified with the input quantity FNY(I). The load NY(I,L) 
in plate element I for load case L is given by 
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NY(I,L) = FNY(I) • NY(L) 


in which NY(L) is the applied transverse load for load case L. 

If PNY(I) 0, NY ( C; js calculated accoruing to modeling 

rules. Therefore, if a zero value of NY (I) is desired for plate 

element I and if the modeling rules give a nonzero value, 

enter FNY(I) as a very small number. 

Shear load, .- The modeling rules in PASCO automatically 

^i 

specify that a singly-connected substructure carries no shear 3oad. 
If two or more doubly-connected substructures join the same two 
nodes, the shear load is distributed between these doubly- 
connected substructures according to the stiffnesses of the shear 
paths. Details of the stress analysis are presented in reference 1. 

Just as in the case of transverse loads, the load NXY(I,L) 
in plate element I fer load case L can be specified by 

NXY(I,L) = FNXY(I) • NXY(L) 

where NXY (L) is the applied shear load for load case L. If 
FNXY(I) = 0, NXY(I) is calculated according to modeling rules. 

PLOTTING 

Plotting is controlled by the parameters AMP, IP, and PAGE 
which are defined in a subsequent section denoted Input. 

When plots are obtained, special attention should be given 
to the first two plots to verify the correctness of the model. 

The first plot shows the repeating element with individual layer 
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thicknesses and element offsets. The second plot shows the 
geometry of each plate centerline for the complete structure. The 
first and second plots are generated independently. The second 
plot/ rather than the first plot, represents the geometry 
actually analyzed by VIPASA. 

The smeared orthotropic stiffnesses and 

calculated within PASCO using the same information as that used 
to gerierate the plot of the repeating element. For that reason, 
if and s.se to be calculated correctly and if subsequent 

analyses based on and are to be correct, it is 

generally necessary for the repeating element shown on the first 
plot to be correct. Analyses that use and include 

analyses involving an loading with SHEAR 0 , a lateral 

pressure loading, a bow- type initial imperfection, temperature 
with ITHERM = 1, and sizing problems involving upper or lower 
limits on ^j^l* 

If plate numbers are reused or if the same number is used for 
mors than one plate, some portion of the repeating element may be 
missing. In the case of an open section, such as a blade- 
stiffened panel, if the Z coordinate is correct for each plate 
element shown in the repeating element and if the complete 
structure shown in the second plot is correct, then and 

are calculated correctly and all analyses are correct. In the 
case of a closed section, such as a hat-stiffened panel, all 
plate elements making up the closed section must be shown 
correctly for the stiffness to be calculated correctly. 
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Analyses that do not use ^11 or D ^2 do not require a correct 
repeating element plot. 

The circular symbols in the second plot indicate node 
locations in the final structure. Mode shapes, which are shown in 
subsequent plots, are interpolated from displacement and rotations 
at these nodes. For a plot to show complex and detailed buckling 
behavior, it may be necessary to introduce a large number of nodes 
through ICARD input. It is emphasized, however, that additional 
nodes do not increase the accuracy of any PASCO eigenvalue analysis. 

For shear and/or anisotropy, two mode shapes are given for 
each eigenvalue corresponding to f^ and f^ in equation (37) , 
reference 1. 

The final plot shows FACTOR as a function of half-wave- 
length X . 

COMPUTER PROGRAM DEFINITIONS AND ORGANIZATION 

In this section, the computer program is described in terms 
of input, output, and program organization. 

Input 

PASCO input is in namelist form with one title card. The 
order of the input is as follows: 

Title Card 
Namelist CONDAT 
Namelist PANEL 
Namelist MATER 
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Input that is considered to be essential is identified by an 
asterisk to the left of the FORTRAN name. 

•V 

Namelist CONDAT .- Namelist CONDAT contains many of the input 
parameters nomally used by the optimization code CONMIN (refs. 5 
and 6) . Parameters that normally appear in CONMIN input data 
^ -t that are not included in namelist CONDAT are X, VLB, VUB, 

ISC, NFDG, FDCH, FDCHM, NDV, NCON, NSIDE, ICNDIR, CT, CTMIN, 

THETA and N1 to N5. Most of these quantities are either 
generated internally based on other input data or are fixed 
because of the approximate analysis approach. For example, 

ISC(N2) is set equal to 1 internally because the Taylor series 
expansion of the constraints used by CONMIN is always linear in 
the sizing variables. 

All of the CONDAT variables which can be input have default 
values that handle normal problems, so it is only necessary to 
input the namelist card. Variables that can be input are defined 
as follows: 

ABOBJl 

ALPHAX 


Fractional change in the objective function 
attempted during initial step in first 
one -dimensional search. May change during 
optimization. Used for interpolation. 
Default is 0.1. 

Maximum fractional change in any X(I) 
attempted during initial step in first 
one-dimensional search. May change during 
optimization. Used for interpolation. 
Default is 0.1. 
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CTL 


CTLMIN 


DABFUN 


DELFUN 


IPRINT 

0 ; 

1 ; 

2 ; 

3-5; 

ITMAX 

ITRM 


Initial constraint tolerance to identify 
active constraints. IF G(J) ^ CTL the 
constraint is defined as active. CTL 
is a negative number. During an optimiza- 
tion the absolute value of CTL is reduced. 
The minimum value of the absolute value of 
CTL is CTLMIN. Default for CTL is -0.1. 

Constraint tolerance to identify violated 
constraints. If G(J) > CTLMIN the con- 
straint is defined as violated. CTLMIN 
is a positive number. Default is 0.000001. 

Convergence tolerance on absolute change 
in the objective function (OBJ) . If in 
ITRM consecutive iterations ABS(OBJfJ) - 
OBJ(J-l)) < DABFUN and no constraints are 
violated, the minimization process is 
terminated. If constraints are violated, 
five iterations are required to terminate. 
Default is 0.000001. (This is a convergence 
test on a "ingle sizing cycle, not a 
convergence test on PASCO.) 

Convergence tolerance. Same as DABFUN 
except comparison is on relative change in 
objective function ABSd.O - OBJ(J-l)/ 
OBJ(J)) instead of absolute change. 

Default is 0.00001. 


Print control for CONMIN. 
is 1. 


Default value 


Print nothing. 

Pring initial and final function information. 

Print all of the above plus control 
parameters. Print function values and X- 
vector at each iteration within CONMIN. 

Increasing amounts of print, used primarily 
for debugging. 

Maximum number of optimization iterations 
performed by CONMIN during each sizing 
cycle. Default value is 15. 

Number of consecutive iterations to 
indicate convergence by DABFUN or DELFUN. 
Default is 8. 
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Linear objective function identifier. 

If the objective function, OBJ, is known 
to be a strictly linear function of the 
sizing variables, X(I), set LINOBJ = 1. 

If OBJ is not a linear function of X(I) , 
use the default value, o. 

Scaling control parameter. 

NSCAL.EQ.O; Do not scale the sizing 

variables . 

NSCAL.GT.O; Scale the sizing variable 

every NSCAL iterations. 
Variables are normalized 
so that scaled X(I) = X(I)/ 
ABS (X(I) ) . 

When using the scaling 
option, it is desirable 
that NSCAL = NDV+1 where 
NDV is the number of 
sizing variables. Default 
is 0. 

Namelist PANEL .- Namelist PANEL provides information 
concerning the initial panel configuration, design conditions, 
sizing variables, and boundary conditions. The integer sub- 
scripts have the following meanings: 


LINOBJ 


NSCAL 


SUBSCRIPT 

I 

II 

IW 

J 

JJ 

K 

L 

M 

MA 

N 


REFERS TO 

plate I 

sizing variable II 
IWALL number 
general ply nunUjer J 
general integer array 

layer in a given wall 
load case L 

eigenvalue wave number in longi- 
tudinal direction 

material number MA 

equation number N that links the 
sizing variables 
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The program is dimensioned to allow large enough numbers of input 


quantities for most practical problems. If there is need, 
dimensions can be increased as indicated in a subsequent section 


entitled Program Operation. After the data are read and 
processed, it is compressed into a blank common block so that 
core size is only that required for the current problem. All 


variables are defaulted to zero unless otherwise specified. 

AB(I,N), AT(J,N) Matrices which can be used to define a 
ATHET(J,N), AC(N) general linear relation between the 

variables as follows AB(I,N) • B(I) + 

AT(J,N) • T(J) + ATHET(J,N) • THET(J) = 
AC(N). These linking equations are used 
to calculate the values of B's, T's, 
and THET's that are considered to be 
dependent variables. Dependent variables 
are identified by using input values = 1.E30. 


AMP Plotting parameter that gives the amplitude, 

in inches, of the buckling mode shape. 
Default, AMP = 0.5 times PAGE. 


AllL(L), AllU(L) Lower and upper bounds, respectively, 

prescribed for smeared orthotropic 
extensional stiffness of panel. 


A33L (L) , A33U (L) 
* B(I) 

BL(I) , BU(I) 


Lower and upper bounds, respectively, 
prescribed for smeared orthotropic shear 
stiffness of panel. 

Width of plate element I. If input as 
negative, will be used as a sizing variable 
with positive value. 

Lower and upper bounds on B(I). Default 
value of BU(I) is I.EIO. 


BMOV(I) 


Governs move limits for sizing variables 
associated with plate widths. Move limits 
for each sizing cycle are generated 
internally and are given by 

VLB(II) = X(II) - DVMOV(II) • (SFACTR)"”^ • X(II)^^^^ 
VUB(II) = X(II) + DVMOV(II) • (SFACTR)”'’^ • X(II)^j^^^ 
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where : 


• VLB (II) and VUB(II) are the sizing \ 

variable lower and upper bounds used by 

CONMIN in a sizing cycle. 

• X(IT) are the values of the sizing 
variables at the beginning oI a sizing 
cycle. 

• DVMOV(II) is equal to the appropriate 
corresponding element in one of the three 
vectors BMOV(I) , TMOV(J) , or THMOV(J) . 

For example, if X(3) = B(5), then DVMOVO) = 

BMOV(5). If BMOV(I), TMOV(J) , or THMOV(J) 
are not read in, the input scalar SMOVE i» 
used for DVMOV(II) . 

• SFACTR is an input scalar 

• n is the sizixig cycle number 

• X(II)j^nit the initial (input) values 
of X(II) . 

Each move limit involves a change that is a 
percentage of the initial value of the 
associated sizing variable. One of the 
objectives of the move limit strategy is 
to reduce the move limits as the sizing 
progresses. Overall lower and upper 
bounds on the sizing variables, defined 
by 'L(I), TL(J), THETL(J) and BU(I), 

TU(J), THETU(JK respectively, override 
the lower and upper bounds VLB and VUB 
for d sizing cycle. 

CLAM(M) Parameter used with buckling or frequency 

constraints to specify margin of safety 
for each half-wavelength A = EL/M. 

Default = 1. 

G = 1 . F(M) 

CLAM(M) 

If the input parcimeter SHEAR is not equal 
to zero, orthotropic plate buckling loads 
are calculated and applied as a correction 
to the M=1 VIPASA buckling load. The adjust- 
ed shear analysis and the df;finition of F 
are discussed in referenc ■ 1 . 
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CONVl 


CONV2 


DELDV 

DELT, 


DELM 


Eigenvalue convergence criterion used 
during eigenvalue analysis to determine 
important buckling half-wavelengths. 
Eigenvalues are calculated to an accuracy 
of one part in CONVl. Default CONVl = 500. 

Eigenvalue convergence criterion. Once 
the important buckling wavelengths are 
identified using CONVl as a relatively 
rough convergence criterion, eigenvalues 
that are to be used as constraints are 
calculated to an accuracy of one part in 
CONV2. Derivatives of these eigenvalues 
are also calculated with CONV2. 

Default = 500,000. 

DELS Parameters used in determining incremental 

DELTH value of sizing variables for numerical 

derivatives . 

AX = DELDV • X(II). .^ + DELX 

init 

where 

are the initial (input) values of 

X(II) . 

DELX = DELB if X is an element width, 

DELX = DEI.r if X is a ply thickness, and 

DELX = DELTH if X is a ply orientation 
angle . 

DELDV, DELB, DELT, and DELTH are input 
quantities . 

DELX is assigned the proper value internal- 
ly. Default, DELDV = 0.0005. 

Controls range of eigenvalues that are 
examined for numerical derivatives. To 
speed convergence of the perturbed eigen- 
value solution, the range of F examined 
for the perturbed solution is 

F(JJ) t AF(JJ) 

where 

AF(JJ) = DELM • DELDV • FACTOR (JJ) 
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DllL(L) , 

ECC(L) 

* EL 

FNXY (I) , 


* FREQ(L) 
FSTIFF 

GRANGE 


and F(JJ) i' the solution for the unper- 
turbed values of the sizing variables. 

(See fig. 28, ref. 1) Tt should be noted 
that AF varies with the finite difference 
increment via DELDV. However the finite 
difference increment also depends upon 
DELB, DELT, and DELTH. 

Increases in DELB, DELT, or DELTH may 
require increases in DELM. Default, 

DELM = 10. 

DllU(L) Lower and upper bounds, respectively, 

prescribed for smeared orthotropic bending 
stiffness of panel. Stiffness resists 
panel bending moment M . 

Amplitude of bow-type imperfection at 
panel midlength for load case L. 

Panel length L. 

FNY(I) Optional input that can be used to override 
the panel internal loads N^y and Ny 
that are calculated automatically by the 
program according to modeling rules. The 
value of Nxyj^ or Nyj^ in any plate i 
is given by 

N (L) = FNXY (I) • NXY(L) 

XYi 

N„ (L) = FNY(I) • NY(L) 

J i 

where NXY and NY are the applied (input) 
loads for load case L. If a user wishes 
to specify a zero value of Njjyj^ or Ny . 

in an element, when the modeling rules give 
a nonzero value, he should input a very 
small number (l.E - 10) rather than zero. 

Frequency requirement for load case L 
in cycles per unit of time. 

Parameter used to specify half-wavelength 
A used by VIPASA to evaluate stiffnesses 
^22» ^33' ^ 22 ' ^ ~ f’STIFF • EL. 

Default, FSTIFF * 10. 

Constraint deletion parameter. Indicates 
range of constraint values re+ained for 
Taylor series expansion and CuNMIN. For 
lower bound constraints, the minimum 
value of G retained is 
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G 


1 


GRANGE 


For upper bound constraints, the minimum 
value of G retained is 


* HOARD (JJ) 


I BAY 


* IBC 


ICARD(JJ) 


^ ~ GRANGE " ^ 

Default, GRANGE = 2. 

Integer array used to create substructures 
and tc define geometry of repeating element. 
The HOARD input is used to rotate plates 
and substructures, offset plates and 
substructures, and join plates and 
substructures. A complete discussion of 
HOARD input is presented in th«^ section 
entitled Modeling. 

For repetitive structures described by 
ICARD input, IBAY is an integer which, 
when added to any node number, produces the 
correct corresponding node number in the 
next repetitive structure. Used only when 
lOREP > 1. Default is nn-1, where nn is 
the maximum node number in the repeating 
element described by the ICARD input. 

•nn - 1 is the correct value for IBAY 
if r< peating elements are connected only 
at onn node.) 

Controls boundary conditions along edges 
y = constant. Not used when ICARD input 
is used. See page 7. 

= 0 free edge 

= 1 simple support 

= 2 symmetry 

= 3 clamped 

Optional integer array used to assemble 
final structure from plates and sub- 
structures created by HOARD. Used when 
more detail is desired in mode shapes or 
when repeating element cannot be described 
by HOARD. A complete discussion of ICARD 
input is presented in the section entitled 
Modeling. 
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ICREP 


IJB 


lORTH 


IP 


I PAST 


TTHERM(L) 


Can be used with ICARD input when structure 
is repetitive. The structure defined by 
ICARD will be repeated ICREP tiroes. 

Controls diagnostic printing. If set equal 
to 1, inforroation concerning negative 
roots is printed during every VIPASA 
iteration. 

If set equal to 0 , anisotropic bending 
stiffnesses are ignored. If set equal to 1, 
anisotropic bending stiffnesses are 
retained. Default = 0 when Nj^y = 0; 

*1 when Nyy / 0. 

Plot and print indicator. 

= 0 no plots 

= ±1 plot of undeformed structure. This 

value of IP is used to obtain only 
a stress and stiffness analysis. No 
buckling analysis is carried out 

= ±2 plot of undeforroed structure plus 

plot of buckling roodes superiroposed 
on undeforroed structure 

= ±3 plot of buckling mode only 

If negative sign is used, the output 
includes coordinates of each plate eleroent 
and substructure defined by HCARD and the 
coordinates of every plate appearing in the 
assembled structure. 

Maximum number of iterations allowed in 
VIPASA to determine eigenvalue for one- 
half wavelength. Default = 100. 

Specifies treatment of bending moment 
caused by temperature and transverse load. 

If ITHERM^l, the panel is allowed to take 
on a bow, the magnitude of which is 
calculated to produce zero bending moment 
in a panel loaded only by temperature and 
transverse load. If ITHERM»0, the panel 
remains flat, and the bending moment 
caused by temperature and transverse load 
is retained. See also discussion of 
bending loads in reference 1 and examples 
in reference 7. 
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IWALL(I) 


* 


JDER 


JPRINT 


* KWALL(K 


LINK 


Wall nvunber for plate I. See section \ 

entitled KWALL and IWALL Input. Default, 

IWALL(I) = I. 

Derivatives of buckling constraints can be 
calculated with a rapid, one-iteration 
approximation or a complete analysis. 

JDER is an integer which indicates the 
number of the sizing cycle at which the 
derivative calculation changes from 
approximate to complete. For example, 
if JDER = 5 and MAXJJJ = «J, the one- 
iteration approach would be used for 
sizing cycles 1 to 4, and complete buckling 
analyses would be used for sizing 
cycles 5 to 8. Default, JDER = MAXJJJ + 1, 
which means that approximate derivatives 
are used for all sizing cycles unless 
JDER is specified otherwise. 

Print indicator. The higher JPRINT, the 
greater the amount of printout. 

1 buckling solution printed only for 
initial design and final design 

2 buckling solution printed for every 
sizing cycle 

3 all buckling solutions printed, 
including those used to determine 
derivatives 

IW) A sequence of integers indicating the 

stacking sequence for wall number iW. 

See section entitled KWALL and IWALL 
Input. 

Indicates whether automatic linking is to 
be used for B(I). If LINK = 0, automatic 
linking is not used. If 

(1) LINK >= 1 and 

(2) plate elements I and J have the 
same wall construction and 

(3) B(J) « B(I) 
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\ 


MAT(J,L) 


* MAXJJJ 


MAXL 


* MINLAM 


MSEL 


* MX(L) 


then automatic linking causes B(J) to 
equal B(I) throughout the sizing. A 
complete discussion of linking is presented 
in the section entitled Linking. The 
default value of LINK is 1. 

Specifies the material number used in the 
jth ply for load case L, This parameter 
is a function of load care to allow material 
properties to be different for different 
load cases. Default, MAT(J,1) = 1, 
and MAT(J,L) = MAT(J,1) for L > 1. 

Number of sizing cycles used in an overall 
synthesis. No convergence criteria are 
provided to reduce the number of sizing 
cycles. When MAXJJJ = 0, the code carries 
out only an analysis of the initial 
design . 

Maximum number of values of half-wavelength 
A for which buckling or frequency 
constraints are calculated. See 
reference 1, section entitled Identifying 
critical buckling and frequency constraints. 
Default, MAXL = 15. 

Specifies smallest half-wave length for 
which buckling loads are examined, 

EL/MINLAM. Default, MINLAM = 1. 

Selects VIPASA mode shape to be plotted. 

For certain problems, numerical ill- 
ccnditioning can occur in calculation of 
mode shapes. An attempt is made to remove 
this ill-conditioning by setting certain 
nodal displacements equal to zero. 

MSEL B 1 specifies mode shape with 
ill-conditioning removed. MSEL * 2 
specifies mode shape as originally 
calculated. Default » 1. 

Applied bending moment per unit width for 
load case L. Should be used when loading 
is primarily bending and it is desired to 
have bending forces included in the 
calculation of the A « EL buckling 
mode. Bending forces caused by M^ are 
not influenced by inplane loads. See 
equation (20), reference 1. 
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NEIG(M) 

NLAM(JJ) 


* NOBAY 
NOSUB 

* NX{L) , NY(L) , 
KXY(L) 

PAGE 

* PRESS (L) 
SFACTR 

SHEAR 


SHOVE 


Number of eigenvalues determined at half- 
wavelength A = EL/M. Default, NEIG = 1. 

Specified values of M (A = EL/M) for 
which buckling or frequency eigenvalues 
are calculated. The program automatically 
calculates the eigenvalue for M = 1. 

The pr'^grcim also calculates eigenvalues 
for other values of M based on logic 
discussed in reference 1, section entitled 
Identifying critijal buckling and 
frequency constraints. With NL^M input, 
the user can specify additional values of 
M for which eigenvalues are calculated. 

Number of major substructures in final 
structure. Default, NOEAY = 1. 


Number of repeating elements in a major 

substructure is equal to 
Default, NOSUB = 1. 


Applied inplane longitudinal, transverse, 
and shear load per unit width of panel 
for load case L. 


Controls size of plot. Default = 1. 

Uniform lateral pressure loading for 
load case L. 


Governs the rate at which move limits are 
decreased. See definition of BMOV(I). 
Default value of SFACTR is 0.8. 

SHEAR ^ 0 causes orthotropic plate buckling 
loads to be calculated and applied as a 
correction to the M * 1 VIPASA buckling 
load. The value of the twisting stiffness 
used in the calculation of the orthotropic 
plate buckling loads is the product of 
SHEAR and the calculated value of the 
twisting s+^iffness. See reference 1, 
section er .itled Adjusted Analysis for 
Shear Buckling. See also reference 8. 


If BMOV(I), THMOV(J), or TMOV(J) are not 
read in, the scalar SMOVE is used in their 
place. See discussion of BMOV(I), Default 
value of SMOVE is C.2. 
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T(J) 


* TEM(J,L) 

* THET(J) 


4" Vi 

Thickness of J ply. If input as negative, 
will be used as a sizing variable with a 
positive value. 

Change in temperature of ply j for load 
case L. 

Fiber orientation angle of ply. If 
input as negative will be used as a 
sizing variable with a positive value. 


THETL(J), THETU(J) Lower and upper bounds, respectively, on 

THET(J). Default value of THETO(J) 
is I.EIO. 


THMOV(J), TMOV(J) Governs move limits for sizing variabl-=;s 

associated with fiber orientation angle and 
ply thickness. See definition of 
BMOV(I) . 

TL(J) , TU(J) Lower and upper bounds, respectively, 

on T(J). Default value of TU(J) is l.ElO. 


Namelist MATER .- Namelist MATER provides material properties, 
indicators for material strength criteria, and allowables used in 
material strength criteria. Subscript MA refers to material 
numbej . Any orthotropic material may be used. Special emphasis 
has been placed on composite materials with direction 1 as the 
fiber direction and direction 2 as transverse to the fiber 


direction. 


El (MA) ,E2 (MA) Young's modulus of orthotropic material 

in 1 and 2 directions, respectively. 

EJ 2 (MA) Inplane shear modulus of orthotropic 

material in 1-2 coordinate system. 


* ANUl (MA) 

* RHO (MA) 


Poisson's raiiio of orthotropic material 
defined by ANUl * E2 = ANU2 • El. 

Mass density of material. 


* ALFAl (MA) 


Coefficitr.t of thermal expansion of 
orthotropic material in 1 direction. 
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ALFA 2 (MA) Coefficient of thermal expansion of 

orthotropic material in 2 direction. 

ALLOW (JJ,MA) Constants used to define material failure 

criterion and allowables used in criterion. 
ALL 0 W( 1 ,MA) is used to select failure 
criterion. A value of 1 indicates an 
allowable stress criterion; 2 indicates 
allowable mechanical strain; 3 indicates 
the Tsai-Wu criterion. ALL 0 W( 2 ,MA)- 
ALL 0 W( 7 ,MA) are input values of the 
allowables defined in the following table. 


TABLE 2 .- INPUT VALUES FOR ALLOW (JJ,MA) 


JJ 

ALLOW (JJ,MA) 

1 

1 Stress 

2 Strain 

3 Tsai-Wu 

2 

comp. 

(£3 - Uj^AT) comp. 

0-^ comp. 

3 

a, ten. 

X 

(£3 - a^AT) ten. 

O3 ten. 

4 

02 comp. 

(£2 - 02 AT) comp. 

0 2 comp . 

5 

O2 ten. 

(£2 ~ 02^"^) ten. 

02 ten. 

6 

■^12 

^12 

■^12 

7 



^12 


In the table, a 21 ^nd allowable stresses; 

E2» and allowable strains; and U2 are coefficients 

of thermal expansions; AT is the change in temperature; ^12 is 
a quantity that appears in the Tsai-Wu material strength criterion 
(see eq. ( 49 ), ref. 1); and ten. and comp, mean tension and 
compression. The allowable stresses and strains are entered as 
positive in compression and negative in tension. 
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Output 


PASCO output begins with an organized listing of input data 
and default values that define the problem to be worked. There 
follows a description of each laminate wall defined by the input 
KWALL. The description includes the materials, thicknesses, 
orientation angles, mass, and laminate elastic properties. Then 
the. overall panel properties, including panel geometry, load 
distribution, mass, and smeared orthotropic stiffnesses, are 
presented. Results of the buckling and stress analyses are 
printed in various levels of detail depending upon the print 
indicators. The analysis of the initial design is then summarized 
in a table which gives the values cf the critical constraints and, 
if struc ;ural sizing is desired, the derivatives of these con- 
straints with respect to each active sizing variable. This table 
contains all data used to generate the Taylor series expansions 
of the constraints. 

The constraint table of the initial design is followed by 
CONMIN output describing the first sizing cycle. The panel 
configuration obtained by CONMIN is then analyzed, and a new 
constraint table is printed. The print pattern of complete 
analysis and CONMIN sizing is repeated MAXJJJ times, then the 
final design is defined in detail. 

Several exeunples of PASCO output are presented in the 
section entitled EXAMPLES. 



Progrcun Operation 


The PASCO computer program was developed on a CDC computer. 
Interim versions have been made operational on IBM computers. 

From this experience, an atLempt has been made to minimize 
changes required to make the progreun operational on IBM and 
other computers. In many cases, double precision is required 
on IBM computers or other computers with similar-length words. 

Organization . - The general program organization is outlined 
in figure 14 . There are a number of comment cards distributed 
throughout the soiirce code that provide additional details 
regarding program organization and logic. 

Dimensions The program is dimensioned to allow large 
enough numbers of input quantities for most practical problems. 

If there is need, dimensions can be increased (or decreased) by 
changing common block BL5 and associated variables that appear in 
PROGRAM COMMON (CDC overlay version of PASCO) . Instructions for 
redimensioning are given in detail in the comment cards in 
PROGRAM COMMON. Variables that can be redimensioned are: 

• Maximum ntimber of plate elements , NOBZ 

• Maximum number of layers in a given wall, IKMAXZ 

• Maximum number of unique layers throughout the structure, 
MAXNLZ 

• Maximum number of values of X for which eigenvalues are 
calculated, MNLAMZ 

• Maximum number of load cases, LOADZ 
56 



} 

I- 


( START ^ 


I 


CAa OVERIAY 3 

READ AND PROCESS INPUT DATA 

X 


CAa OVERLAY 4 
PRINT INPUT DESIGN 



CAa OVERIAY 5 
SETUP VIPASA INPUT ON TAPEl 


CAa OVERLAY 1 
PERFORM VIPASA ANALYSIS AND OTHER ANALYSES. 
INCLUDING DERIVATIVES 


CAa OVERIAY 2 
OPTIMIZE WITHIN MOVE LIMITS 



CAa OVERLAY 2 
PRINT CONSTRAINT ARRAY 



Figure 14.- PASCO program flow diagram. 
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• Maximxim mimber of integers appearing in HCARD and ICARD 
input, JJHZ and JJIZ, respectively 

® Maximum number of materials , IDUMJZ 

• Maximum number of linking equations, MIZ 

• Maximum number of constants used to define material 
strength criteria, NALLOW 

Diagnostics . - An attempt has been made to provide diagnostics 
for common input errors that lead to inconsistent or incomplete 
problem definition. Undoubtedly some errors will not be detected 
by the diagnostics. Also, unreasonable proportions can lead to 
numerical ill-conditioning that causes program dumps. 

Plotting . - Because individual computer centers generally 
have their own plotting commands, the plotting command subroutines 
that are used in the version of PASCO operational at Langley 
Research Center are not included in the version of PASCO that is 
made available to other users. The subroutines that execute the 
plotting in the Langley version of PASCO are described with 
comment cards in the source listing in sufficient detail that 
user subroutines Cin be substituted to perform the same functions. 
Additional information on plotting subroutines is given in 
reference 4. 
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EXAMPLES 


Several examples are presented to illustrate in detail the 
features of and modeling techniques used in PASCO. For each 
example, the following information is presented: (1) a computer 

drawing of the repeating element of the initial design, (2) a 
computer drawing of the cross section of the entire panel with 
nodes, (3) a figure showing the design loading, (4) a listing of 
the input, and (5) selected program output. Taken as a whole, 
the selected output for each exeunple encompasses the major portion 
of the output of the program. The numbers on each repeating 
element are the original plate element numbers. Units are not 
used in input or output quantities in the computer program and are 
not shown on the figures. 

Example 1, Unstiffened Plate 

The first example is a laminated, unstiffened plate loaded 

in longitudinal compression as shown in figure 15. The stacking 

sequence shown in figure 15(c) is created with the input T, THET, 

and KWALL presented in figure 16. The stacking sequence is 

(+450, -45°, 90°, 0°) with corresponding thicknesses 

s 

(.02, .02, .02, .01)g. The negative signs on the T and THET 
input data indicate that the sizing variables are T(l), T(2), 
T(3), THET(2), and THET (3) . For this simple example which has 
only one plate element, the default value IWALL(1)=1 is correct. 
HCARD input is not required because the entire structure is one 
plate element. 
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(b) Entire panel, with nodes 


///// 



(c) Panel with design loading 

Figure 15.- Exeunple 1, unstiffened laminated plate designed 
for longitudinal compression. 
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UNSTIFFFNED PLATE 

♦CONOAT 

IPPINTs?, 

% 

«PANEL 

MAXJJJsl* 

B(l)sP.» 

-.0?* -,0P* 

THFT(n*0.» -AS., -90.* 

KWALLn»l)«2.-2* 3. 1* 

Nxn)sioop,. 

IRC*l* 

EL»30., 

IPs?. 

H1NLAH=30. 

% 

umatfp 

Ein)si9.F6* F?(1 )s1,89E6. E1?( 1 ) s, 93E6. RHO ( 1 ) s.05Tl , 

AlFAl n )=-,00*)E-6. ALFA2(1 )s21,8E-6. 

ALLOW(1«1)s2, .004* -.004* .004* -.004* .01* 

% 


Figure 16.- Input for example 1, unstiffened laminated plate 
subjected to longitudinal compression. 
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The output selected for this example is simply output related 
to the model. It describes the wall construction and gives the 
wall properties for the one wall used in this exaunple. The output 
is presented in figure 17. In other analysis and sizing problems, 
a similar set of wall property information is presented for each 
unique wall. In PASCO, where all walls are balanced and 
symmetric, ZREF is half the total thickness. The quantities 
E12, denoted LAMINATE PROPERTIED, are given by 

^11 ~ ^^ 12 ^ ^^22 
Total Thickness 


El , £2 , and 

El = 


E2 = 


^22 ~ ^^ 12 ^ ^^11 
Total Thickness 


E12 = 


'33 


Total Thickness 



*»*IL 


MIIMRCP 


1 


MAT 

T 

TMCT 

ten 

1 

.0200000 

OS. 00000 

O.OOOOO 

1 

.0200000 

•05.00000 

O.OOOOo 

1 

.0200000 

90.00000 

o.eoooo 

1 

.0100000 

0.00000 

O.OOOOo 

1 

.0100000 

0.00000 

O.OOOOo 

1 

.0200000 

90.00000 

O.OOOOO 

1 

.0200000 

-4S.00000 

O.OOOOo 

1 

.0200000 

45.00000 

O.OOOOO 

ZftCF FOP WALL 

NO. !• 

.7000E-01 



MSS PCR UNIT ARE* Of MALL NO. 1> .T994E-0? 


A-MATRIA 

FOR MALL 

NO. 1 

DmMATPIX fop wall 

NO. 1 

CXT 


.97992*0* 

.40**E*0* 

.03S2E-OS 

•1415C#04 •9079F«O3 

•3499e«03 

Os 


.40**F*0» 

.1325E*07 

.2650E.0* 

•9«79F«03 •1703C«04 

•3495C^03 

Os 


.•392C-0N 

.2«50E>0* 

.4540E*0« 

,34$5e«03 •34S9F«03 

•1049r^04 

Os 





LANINATE RRORERTIES 



fl 

C2 

C12 

NUl 

NU2 

ALFAl 

«LF*2 

610«31E*07 

• S?6206F407 ,324264C*07 .30**^674 . 

4149314 . 

35247aC«0B 

.1S737«E 


Figure 17.- Sample output for example 1, unstiffened 

laminated plate designed for longitudinal compression. 
Output gives wall property information. 
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Example 2 , Blade-Stiffened Panel 

The second example is a blade-stiffened panel loaded in 
longitudinal compression as shown in figure 18. 

The input data T, THET, and KWALL shown in figure 19 are 

used to create three walls. For example, wall 3 has a stacking 

sequence of (+45°, -45°, -45°, +45°, 0°) with corresponding 

s 

thicknesses (.004, .004, .004, .004, .04) . The input IWALL(I) 
assigns the wall number of plate element I. For example, 
plate element 4 has wall 2. 

The panel geometry is modeled with one line of HCARD data 
consisting of three data sequences. The input NOSUB=3 produces 
four repeating elements per major substructure or bay and 
determines the spacing of the nodes. The input NOBAY=4 produces 
four of these bays for the complete panel. The panel, therefore, 
has sixteen stiffeners. 

The negative signs on B(l) and B(2) cause the stiffener 
spacing to be variable during panel sizing. Since the number of 
stiffeners is fixed, the panel width becomes variable for this 
example. (Linking equations could have been used to fix B(l) + 
B(2).) Automatic linking is used to cause B(4) to equal B(2) 
and B(5) to equal B(l) during sizing. 

The output chosen for this example shows a portion of the 
organized listing of input data and default values that define 
the problem. This output, which is presented in figure 20, 
shows the plate element widths B(I), the layer thicknesses T(J), 
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and the layer orientation angles THET(J), together with their 
lower and upper bounds for sizing. 


I 



I ? 4 5 



(a) Rspeating elernent 


4 


{f)) Entire panel, wiln nodes 



^ ' / / >' ^ / y y / 

• ^Uix' 

iO Mth deshjfi (cadin<^ 

Figure JH.- Example 2, Ijiudo-stif fened panel with sixteen 
stiiteriorri, designed for longitudinal compression. 
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*«•*«**•«• BLAOF-STIFFENED PANEL ♦•***#♦*♦#♦••*•*• 

fCONOAT 

$ 

upanel 

HAXJJJs2t 

LINKai, 

Ip5?f 

••2f *l»Sf *?f •ft* 

T(1)3 -.OO^* >.00«* -.003* >•01* ••04* 

THFTU)= 45* 45* 0* 0* 0* 

KWALUI.D* l.-l*-l,l,3* 

KWAU.(I*2) «1*-1*-1,1«4* 

KWAU(1»3) « ?«-2.>7,?«5* 

IWALL(1)»1*2*3«2*1* 

JOFR«l. 

HCARO(l)s 4*-6*3*90,0« 2*121*6. 6*7*lt?«4*»121 *5* 

N0RAY«4* 

NOSURs3* 

IRC*2* 

FL»30,* 

NX(l)s3000.« 

NEIG(1)«4* 

NLAM(1)«1,3»17* 

MINLA»<«30* 

» 

»MATEP 

F1<1)b19.E6* E2(1)>1.R9E6* E12 ( 1 1 «.93E6* ANu 1<1I>.31« RHO ( 1 ) >.nS71 * 

ALEAI (!)«•. 005E-6* ALEA2(l)*31.ftE-ft* 

ALL0W<1*1)« 2* .004* -.004* .004* -.004* .01* 

S 


Figure 19.- Input for example 2 , blade stiffened panel 
subjected to longitudinal compression. 
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Example 3, J-Stiffened Panel 

The third example is the J-stiffened panel shown in 
figure 21. The design loading is longitudinal compression. The 
panel is to be sized to carry that load assuming that the panel 
has either a positive bow or a negative bow. This is an example 
of a two- load-condition problem in which the first load condition 
is Njj = 12000, e = +0.09 and the second load condition is 
= 12000, e = -0.09. These two loading conditions are 
illustrated in figure 21(c). The input variable for e is 
ECC. 

Based on the input presented in figure 22, three walls are 
created. Each of these walls has two materials. For example, 
wall 3 has a stacking sequence of (45°, -45°, 0°, 90°) with 

o 

corresponding thicknesses (0.008, 0.008, 0.016, 0.008) and 

s 

materials (1, 1, 2, 1) . The input quantity IWALL assigns these 
three walls to the fourteen plate elements. 

The repeating element shown in figure 21(a) illustrates 
several modeling features. Offsets created with six HOARD 
sequences provide a more realistic model. The offsets are also 
shown in figure 21(b). The stiffener attachment flange is 
attached to the skin with short plate elements to insure that the 
stiffener attachment flange and the skin directly above the 
flange respond as a unit. For this configuration, both the 
attachment flange and the skin directly above the flange could 
have been divided into two plate elements rather than four 
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(a) Repeating element 



(b) Entire panel, with nodes 



(c) Panel with design loadings 


Figure 21.- Example 3, J-stiffened panel with eight 
stiffeners designed for bow and longitudinal compression 
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J-STTFrCNEO PANEL 
ACONOAT 
S 

»PANEL 

SM0VE»,15f 

MAXJJJ*?. 

LINK«1, 

CLAMs1,3. 

JDFRal, 

NLAM(n» <,10, 
NEIG(1)*4, 

NINLAM*30, 

FL*30, , 

NX(nsl?000«* 12000., 

ECC».09, -.09, 

IP*2, 


THFT(1)>4S.,0.,90,,AS.,0.*90,,0., 

MAT(l,n«l,2,l,l,3,l,?, 

KWALL(l*nEl,.i,2,3, 

FWALL{l,?)aA,-4,5,6, 

KWALL(1,3)>A,.4,7,6, 

I NALL *1 «l,l*?,?,3,3,3,2,ltlt2,2*l, 
HCARO* 0 , — 1 5 , — 4 , 0 , — A , 0 • *6 , 

6,— l6,-S,0,-A,0,-6* 
S,17»2,»15,3,-16» 

6,-18,-12«0,-A,0,~Pt 

0»*19,-13, 0 ,-8, 0 ,-R, 
5,20*10,-18,11,-19, 

A, 21, 8, 7. 6, 

4,-22,21,90,0, 

6 *-23, -22 ,0,0(0, -9, 
2,121,9, 


6 ,- 122 ,- 121 , 10 , 0 , 0 , 0 , 

3,123,23,-122, 

6,15,1,17,20,-123,14, 

N0RAY«4, N0SUB«2, 

ATn,l)>2.,l.,l.,2.,l.,l.,0.,-l., 
AT(1,2)*2.,1.,1.,4,,2.,2.,0.,0.,-1,, 

AR ( 1 , 3 ) >0, , 1 , , 0 , , —1,, 

AT(1,4)«0.,0.,0.,2.,0,,1.,1,,0,,0,,-1., 

A 11L( 11*2200000., A33L (11*450000,, 

Till*-, 04, -.012, -,00A, -.008, -.07, -.008, **•016, 

B(ll*-.5, -.2, .2, 1.E30, .2, -.4, .4, .4, -.4, .2, 
t 

SNA TER 

£1(11*11. E6, E2(’)».«E6, E12(11*.3E6, ANU1(1|*,34, RHO(11*,OS, 
ALFAI (lla.iE-98, ALFA2(1)*1,E-5, 

ALLONd ,1 1*2, .003,-,017,.S3;-.5C5, .01, 

ALLOW(1,21*2..009,-.005..007,-.003,.010, 


3*1. E30, 

2f aPt aSt 


Figure 22.- Input for example 3, J-stiffened panel with eight 
®'tiffeners, designed for bow and longitudinal compression. 
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plate elements.^ Because of the bow, the web of the stiffener 
is modeled with three plate elements to allow the stress in the web 
to vary as the panel bends. (See? fig. 8, ref. 1.) 

When a panel is to be sized to account for a bow, the input 
CLAM(l) should be used to specify a margin on the M=1 buckling 
load. The strategy is to use a combination of CLAM(l) , SM07E , 
and SFACTR that prevents large moves when the constraints are 
highly nonlinear, as they can be when F(l) approaches 1.0. 

(See discussion of bending loads in ref. 1.) It is also advisable 
to use the input NEIG(l) to obtain more than one eigenvalue for 
the M=1 mode. 

This example also illustrates the input for linking equations 
(AT and AB) and stiffness requirements (AllL and A33L) . 

The output selected for this example gives the table of 
linking equations, figure 23(a), and the plate forces and stresses 
for the first load condition , ECC=0.09, figure 23(b). Each 
column in the table of linking equations represents a linking 
equation. The entries above the top row of asterisks are the 
coefficients of the unknowns, the dependent variables. The 
entries between the two rows of asterisks are the coefficients 


^ CAUTION ; Plate elements having narrow widths appear to 
cause numerical instabilities in VIPASA. As a rule of thumb, 
width-to-thickness ratios should be greater than two. If 
unexpected dumps occur, try remodeling with greater widths . 
Plate elements having very long lengths EL can also cause 
numerical instabilities. This problem can usually be solved 
by setting FSTIFF equal to 1.0. 
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vapiablf 


number linking equation no. 



1 

2 

3 

4 

4 

0.0000 

0.0000 

-1.0000 

0,0000 

?2 

-1.0000 

0.0000 

0.0000 

0,0000 

23 

0.0000 

-1.0000 

0.0000 

0,0000 

24 

0.0000 

0.0000 

0.0000 

-1,0000 

2 

0,0000 

0.0000 

1.0000 

0,0000 

IB 

2,0000 

2.0000 

0.0000 

0,0000 

16 

1.0000 

1.0000 

o.oooo 

o.oooo 

17 

1.0000 

1.0000 

0.0000 

0,0000 

18 

2.0000 

4.0000 

0.0000 

2,0000 

19 

1.0000 

2.0000 

0.0000 

o.oooo 

20 

1.0000 

2,0000 

0.0000 

I.OOOO 

21 

0.0000 

0.0000 

0.0000 

1,0000 

AC 

0.0000 

0.0000 

0.0000 

o.oooo 


(a) Table of linking equations 


Figure 23.- Sample output for example 3, J-stiffened panel 
with eight stiffeners, designed for bow and longitudinal 
compression. 
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(b) Plate forces and stresses for ECC-0.09 case 


of the independent variables. The entries in the last row are 
the values of AC. The variables B(I), T(J) , and THET(J) 
are identified by the numbers in the column denoted VARIABLE 
NUMBER. 

It is important to note ho^-e that linking equations are not 
required in a problem involving only analysis — no sizing. For 
analysis applications, input values of B, T, and THET are 
adequate. Linking equations are useful primarily for sizing 
applications or for parameter studies (a series of single 
analyses) where dimensions arc being changed. Assume, for 
example, that an engineer wisiies to carry out a parameter study 
with the J-stiffened panel described in this section and that 
the parameter study involves varying T(l). The offset T(8) of 
plate elements 4, 5, 12, and 13 varies with T(l). The engineer 
can either specify the thickness T(l) and the offset T(8) for 
each case in the parameter study or he can specify T(l) for each 
case and rely on a linking equation that automatically calculates 
the corresponding value of T(8). 

The plate forces in each plate element and the strains and 
stresses in each layer are shown in figure 23(b) . The columns 
denoted PLATE NO give the VIPASA identification for the plate 
elements. Plate numbers 4, 5, and 6 make up the web of the 
stiffener. The effect of the bending moment can be seen by the 
variation of the lead in these plate elements. The strains and 
stresses in each layer are parallel and transverse to the fiber 
direction in that layer. 
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Example 4, I-Stiffened Panel 

The fourth example is an I-stiffened panel loaded by a 
combination of shear and longitudinal compression as shown i.. 
figure 24. 

The modeling for the I-stiffened panel is similar to that 
of the J-stiffened panel discussed previously. In the I- 
stiffened panel example, a fifth linking equation causes the 
stiffener spacing to be equal to 5.0 during panel sizing. 

This example illustrates the adjusted analysis for shear 
buckling in PASCO. The. adjusted shear analysis (rather than the 
standard VIPASA analysis) is selected by setting the input 
parameter SHEAR equal to a number greater than zero. The 
parameter SHEAR is set equal to 1.0 in the input shown in 
figure 25. If the panel were composed of stiffeners with closed 
sections, such as hats, a value of shear less than 1.0 is 
recommended. See definition of SHEAR in section entitled Input 
and in reference 1. 

The program output (fig. 26) selected for this example 
focuses on the adjusted analysis for buckling when the loading 
involves shear. The output consists, primarily, of three VIPASA 
solutions: two solutions obtained with smeared orthotropic 

stiffnesses and one standard VIPASA solution obtained with 
discrete stiffeners. These three solutions are used in the 
adjusted shear analysis and are discussed in the section entitled 
Adjusted Analysis for Shear Buckling in reference 1. All of this 
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N^* 6000 


(c) Panel with design loading 


Figure 24.“ Exeunple 4» I—stiffened p&nel with eight 
stiffeners designed for longitudinal compression and shear 
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Figure 25.- Input for example 4, I-stiffened pemel with eight 
stiffeners, designed for longitudinal compression and shear. 
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(d Stamtoi J VIPASA solution. F. . , and adjusted solution 


NEGATIVF wall numbers indicate anisotropic STT'^fNESS TERMS NEGLECTED 
PLATE FORCES ARE FOR UNIFORM AXIAL STRAIN 


PLATF WALL 


NO, 

NC. 

0REAOTH 

NTV 

NLV 

NSV 

B*NLV 

1 

-1 

1.90000 

0.00 

1612.01 

600*00 

3062.82 

? 

-1 

•300CO 

0.00 

1612.01 

370*99 

483.60 

3 

-? 

•30000 

0.00 

10313.29 

221*01 

3093.99 

4 

••3 

.25000 

0.00 

1751.16 

o*oo 

437.79 

5 

-2 

.40000 

0.00 

10313.29 

o*oo 

4125.31 


(d) Stress •^sultents Ny, and Nxy in each plate element 

Figure 26. -• Concluded. 
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output is for the initial design. Comments added to the output 
use the same notation as reference 1. 

Equation (45), reference 1, is the basic equation used in 
the adjusted shear analysis. That equation is 


F 

d,90 


F 

F s,90 

^s,0 


( 10 ) 


for this example, the quantities that appear in equation (1) are 
given by 


F 

d,90 


.331802647 

.376198540 


.516336069 


= .4554022842 

Since F, „„ is smaller than F on is used for the 

adjusted analysis. 

A multiplier defined to be the smaller of 


F F 

d,90 , s,90 

F 

d,0 d,0 


is then calculated. In this case, the multiplier is equal to 
1.37250947. The adjusted solution is then given by 

Multiplier * FACTOR • [input loading] 
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The multiplier is applied to all eigenvalues for the M=1 mode. 
For this example, the adjusted solution for the lowest M=1 
eigenvalue is 


1.37250947 • .331802832 


6000 (N ) 

600 (N^y) 


which gives ^Xcr ” 2732.41517, - 273.241517. 

The final output data presented for this example are the 
forces in each plate. The shear forces are labeled NSV. Note 
that there are no shear forces in the stiffener. 

See additional shear buckling studies in reference 8. 



Example 5, Hat-Stiffened Panel 

The fifth example is the hat-stiffened panel shown in 
figure 27. The two design loadings are shown in figure 27(c). 

Load condition 1 consists of pure longitudinal compression. 

Load condition 2 consists of longitudinal and transverse 
compression combined with a 200° increase in temperature. 

In the input shown in figure 28, ITHERM(2)=1 allows the panel 
to take on a bow induced by temperature and transverse load. 

This bow is treated as an initial bow. Note that material 2, 
which, in this example, is used for the second load case, has 
lower stiffness values than material 1. Using this approach, 
the effect of temperature on material properties can be accounted 
for. The input also illustrates how the rotation of a plate 
element can be made to depend upon the lengths of various plate 
elements making up the cross section. Such an approach, which is 
based on equation (3) , must be used if the lengths of these 
plate elements change during sizing. 

The output selected for thi.s example is the constraint table 
for the first sizing cycle. Thu constraint table, which is 
presented in figure 29, contains the values of the constraints G 
and the derivatives of the constraints with respect to each sizing 
variable. The first column in the table identifies the constraints. 
For example, the first four constraints are the first four 
eigenvalues for the M=1 mode (X=L) for load condition 1. The 
second column gives the values of the constraints G. In this 
case, the first two constraints are violated, and the second 
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(a) Repeating element 



(W Entire panei, with nodes 


////// 



////// 


////// 

N -4U» 


LOAD CONDITION 1 


////// temperature ■ 200° 

■ 1500 allowed to take 

! ^ ON INITIAL BOW INDUCED 

BY TEMPERATURE AND 
TRANSVERSE LOAD 
LOAD CONDITION 2 


(cl Panel with design loadings 
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'aA. »4- >4. •♦6. .4. »4. (4. 1.F30* »Rr ««1R» 

Ta-,00fi. -.01ft. -.03. ?*1E30. 

THETrAft.O.O. 

AT(l.l)* P.l.O.-l. 

AT(1.?)» ?.0,1.0.-1. 

AB(«5.3)s 1,. 3*0.. -1.. 0.. 6 . . 

MATd.Ds 1. 1. 1. 

MAT(l.?»s ?. ?. ?. 

ITHF.PM(?)=i , 

KWiLLd.Usl.-l.-l.l.?. 

KWALL(1.?J=1.-1. 

KWALL(1.3)*1.-1.3. 

IWALL= 1.7.?.?. 3. 

EL»30. 

NX(U=4S00.. 1500.. 

NV(?)=100.. 

TFM(1.?)=3*?00.. 

:bcs?. 

HCAPns4.1?.?.3.4. 
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4. -14.1?. -11.?. 

4.-15.13.11.?. 

ft. -16. -14. 0.-4. 0.5. 

ft. -17. -15. 0.5. 0.-4. 

6.— IP.- P.0.4.0.4. 

4. 19. 1ft.*;. 17. 

5.70.1 .IP. -19. 10. 

N0BAY=4.N0S0Pa?. 

NFIG(1)=4. 

MINLAP=30. 

% 

SMATEP 

Eldja 71F6. E?d)= 7.39E6. E12d)= .ft^Fft. ANUld)« .314. RHOd'«,0571. 
El(?)s 19E6. F?(?)= ?.00Eft. El2(?)= .ftOEft. ANI11(7)» .314. PHO (?) «.0571 . 
ALFAld)a -.5F-B. A|,FA?(1)= .71AE-4. 

ALEAl(?)s -.5F-«. ALEA?(?)s .?1BE-4. 

ALLOWd.l)*?. .004. -.004. .004. -.004. .01. 

ALL0¥(1.2)*?. .004. -.004. .004. -.004. ,01. 
it 


Figure 28.- Input for example 5, hat-stiffened panel with eight 
stiffeners, designed for longitudinal and transverse 
compression and temperature. 
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two constraints are satisfied. The final seven columns are 
the derivatives of the constraints with respect to the sizing 
variables. The derivatives are given first with respect to B(I) , 
then T(J), then THET(J). In this case, the order of the 
sizing variables is B(l), B(2), B(5), B(ll), T(l), T(2), 

and T(3) corresponding to the negative signs in the input data 
for B and T in figure 28. 



Example 6, Z-Stiffened Panel 

The sixth example is an aluminum Z-stiffened panel loaded in 
longitudinal compression as shown in figure 30. Instead of a 
constraint on buckling, which is used in other examples, in 
this excunple a constraint is placed on vibration frequency. 

Another feature of this example is the method used to model 
the attachment of the stiffener to the plate. Whereas in previous 
excunples the attachment simulates a bond, in this example the 
attachment simulates a seam weld oi a line of closely-spaced 
rivets. 

The input is presented in figure 31. Note that since 
this excunple involves vibration, RHO must be mass density. 

The output selected for this excimple gives the results of 
the eigenvalue analysis for the initial design and the CONMIN 
output for the first sizing cycle* In the eigenvalue output, 
figure 32, FACTOR is applied to both the input load NX=3000 
and the input frequency FREQ=50. For example, in the second 
eigenvalue for LAMBDA=30, FACTOR is equal to 1.84758441. The 
value of NX for this eigenvalue is 1.84758441 • 3000 = 5542.75, 
which is the first entry in the column leibeled NX. The 
vibration frequency for this eigenvalue is 1.84758441 • 50 = 92.38, 
which is not printed out. 

In the CONMIN output, figure 33, the default value 1 is 
used for the print indicator IPkINT. With IPRINT=1, only the 
initial function infonnation and the final results are printed 
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(a) Repeating element 




3 

L — L^X 





o 
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<b) Entire panel, with nodes 





Ny- 3000 


(c) Panel with design loading 


Figure 30.- Exeunple 6 , Z-stiffened panel with eight 

stiffeners, designed for longitudinal compression and 
a minimum frequency constraint. 
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KWALL (1*2)=2, 
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MAT(1*1)»1*1, 
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T«-.l .-.05*l,F30. 
R*“2**2»««»B*,5*«1»?*1»E30* 

PL (3)«.5» 

FPFQ*50.* 

t 

8MATER 

F1*10.5E6, 

F2»10,'jF6# 

E1?»3.98E6* 

ANU1«.32« 

RHOs. 00028138* 

A|.FAlel3,E-6* ALFA?*13.E-6* 
ALLOW*?* .006*-.008*.0J6*-.006*.01 * 
* 


Figure 31.- Input for example 6, Z-3tiffened panel with eight 
stiffeners, designed for longitudinal compression and a 
minimum vibration frequency. 
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FORTRAN OROGRAM FO« 
rONSTRAINFO FUNCTION MTNI^'l 7ATION 


INITIAL function information 
ORJ » ,?4SR55F-0*i 

OFCISIPN VARIABLFR (X-VFCTOR) 

1) ,?OOOOE*ni .ROOOOE^OO .l? 000 F* 0 l ,10000E*00 ,50000F-01 

CONSTRAINT valuer (G-VECTOR) 

1) -.49fll9E»nn -,84758E*00 

final OPTIMISATION information 
ORJ * .?! 1718F-0R 

OFCTSION VARIABLFR (X-VECTOR) 

1) .??OOOE»01 .RO0OOF*O0 ,I080«F*V1 .90000F-0I .'.5000E-0I 

CONSTRAINT VALUER (6-VECTOR) 

I) -.4f069E-0I -.?4795E*00 

THFPF ARf 0 ACTIVE CONSTRAINTS 

TMFRF are 0 violated CONSTRAINTS 

TMFPE are R active SIDE CONSTRAINTS 

OFCIRION VaRIACLRR AT LOWER OR UPPER ROUNOS (MINUS INDICATES LOWER POUND) 
I -? -3 -4 -R 

TERMINATION CRITFRION 

ABS(Opj(n-nRJ( r-I) ) LESS THAN DARFUN FOE 8 ITERATIONS 


NUMBER OF 

ITERATIONS = 




OBJECTIVE 

FUNCTION WAS EVALUATED 

i?. 

TIMES 

CONSTRAINT 

FUNCTIONS WERE 

evaluated 

1? 

TIMES 

gradient of OTjtCTIVF WAS 

calculated 

s 

TIMES 

GRADIENTS 

OF CONSTRAINTS 

WERE CALCULATED 

s 

TIMES 


Figure 33.- Sample output for example 6, Z-stiffened panel. 
Output is from CONMIN program. 
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in CONMIN. In the final optimization information, the quantities 
labeled CONSTRAINT VALUES (G-VECTOR) are approximate values of 
the constraints calculated using the approximate analysis 
approach discussed in reference 1. 



Example 7, Truss-Core Sandwich Panel 

The seventh example is a truss-core sandwich panel loaded 
by a combination of longitudinal compression and lateral pressure. 
The panel and loading are shown in figure 34. The example is 
relatively complex and, in some respects, stretches to the utmost 
the capability of the program in order to illustrate its 
ci--'.neraliti and limitations. 

The input is shown in figure 35. Because of limitations 

associated with HOARD input discussed earlier, ICARD input is 

required to model the panel. Even though HCARD input is not 

adequate to model the repeating element from the point of view of 

correct structural response, HCARD input can be used to model the 

geometry of the repeating element. Because the lateral pressure 

loading causes a bending moment, the geometry of the repeating 

element must be defined by HCARD input in order that D^^ 

and the longitudinal load distribution be correct. 

^i 

The single linking equation states that B(7) = .5 B(l). This 
value of B(7) is used in the rotations specified by the first 
four HCARD sequences. Plate elements 7 to 10 created by these 
first four HCARD sequences are used for ICARD input. The other 
HCARD sequences illustrate one of several ways to construct the 
geometry of the repeating element. 

The output selected for this example is shown in figure 36. 
The upper portion of the figure is the HCARD and ICARD input 
actually used by VIPASA. Note the large amount of ICARD data 



1 



(b) Entire panel, with nodes 


///// 



(c) Panel with design loading 

Figure 34.- Example 7, truss-core sandwich panel with eight 
repeating elements, designed for longitudinal compression 
and lateral pressure. 


96 


I 


ORIGINAL PAGE 
OF POOR QUALITY 
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^=”1 • « “!•* 1«* !•« !•« 1*E30« 
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IWALL=1 t?. ?»?,?. It 
EL=30t 

hCARn=A,-7,?,-7,?t 

4,-P,3t7,2t 

4,-9,4t7t?t 

4,-10t5»-7,2, 

4,-11 tlOt 180,0. 
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3,121,9,-1?, 

4, 13, 7, P, -121 , 

3, 12?, 1, -13, 
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I CAPn=7,l ,3,7,4,1,-909,900, 

7, ?, 3, 9, S, 6, -909, 900, 

5.3.4, P,s,10, 

3.4, -909,900, 

3, <3, -909, 900, 

AB(1,1)=*R, 5*0,, “1,, 

$ 

tMATFR 

El(l)=19.Eft, F?m=l ,P9Efi, E12(1)=.93E6, *^*U1 ( 1 ) *.31 , RHO ( 1 ) », 0S71 , 

AIFA1=-.‘>E-8,ALF42s.?1PF-4, 

ALLOW*?,. 004, -.004, .004 ,-.004 , . 01 , 

< 


Figure 35.- Input for example 7, truss-core sandwich panel 
designed for longitudinal compression adn lateral pressure. 
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here, figure 36, compared with the small amount in the input, 
figure 35. 

At the bottom of figure 36 .'re six smeared orthotropic 
stiffnesses for the panel and quantities denoted ZBAR and AM. 

The stiffnesses are defined in reference 1. The stiffnesses are 
calculated from the repeating element constructed from HCARD 
input, and because of limitations in that model, are not all 
correct. In particular, the values given for the stiffnesses 
A22, A3.1, and D22, which are calculated by VIPASA, are 
meaningless since the repeating element is defined as being 
singly connected. With the present limitations in VIPASA, no 
modeling approach would provide correct values for A22, A33, and 

D22. The values given for All and Dll are correct. The 
value given for D33 is too small. The value is approximately 
the twisting stiffness of the upper triangle of the repeating 
element. Since several pertinent stiffnesses are incorrect, 
the adjusted shear analysis (example 4) should not be used for 
this configuration. It is important to recognize that these 
smeared stiffnesses are not used in VIPASA, and, therefore, 
errors in the stiffnesses do not cause errors in the VIPASA 
buckling analysis. 

The quantity ZBAR is the distance from the reference 
surface to the centroid of the cross section. The reference 
surface is the horizontal plane that passes through the first 
node of the repeating element. The quantity AM is the mass per 
unit area of the panel. Both values are correct. 
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Example 8, Stringer-Stiffened Cylinder 

The eighth example is a composite cylinder stiffened by 
rectangular, integral stringers ("blades") and loaded by 
longitudinal compression. The cylinder and loading are shown 
in figure 37. 

The input is shown in figure 38. HOARD input is used to 
create eight sections of the stiffened cylinder. ICARD input 
is used (and is required) to connect these eight sections to 
form the cylinder. The last three HOARD sequences create the 
repeating element. The repeating element is used to calculate 
the stiffnesses and the value of the uniform strain that is 

used for all plate elements - including the plate elements in 
the cylinder. Note that in order to calculate the stiffnesses 
and e^, the repeating element is flat. Note also that the 
repeating element is not used in the cylinder. 

The input also illustrates that substructure numbers can be 
reused, and that the B(I) used in rotations are the original 
values of B(I) and are not associated in any way with subsequent 
substructures having the Scune value of I. In this example, B(7) 
to B(12) are cosines of the desired angles. As can be seen in 
the input, the skin on the right half of the cylinder is 
constructed of wall 1, and the skin on the left half of the 
cylinder is constructed of wall 3. Wall 3 is identical with 
wall 1 except that the angles are reversed. The cylinder is 
modeled in this way so that after rotation of plate elements 
used on the two halfs of the cylinder, the ply angles in each 



3 



Figure 37.- Example 8, composite stringer-stiffened cylinder 
with radius equal to 40, length equal to 30, designed for 
longitudinal compression. 
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Figure 38.- Input for example 8, stringer-stiffened 
crylinder designed for longitudinal compression. 
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layer are continuous. Differences in the walls (±0 vs. +6) 
would not be detected if anisotropic bending stiffnesses were 
ignored, but with I0RTH=1, as it is in this example, tne 
differences between walls can be detected. 

The output selected for this example consists of data 
describing the plate elements, figure 39(a), and HOARD and ICARD 
input data actually used by VIPASA, figure 39(b). The data 
describing the plate elements gives the VIPASA plate number 
assigned to the plate element widths, and shows at a glance that 
B(3) is the only plate element width that is a sizing variable. 
The HOARD input shown in figure 39 (b) uses the VIPASA plate 
numbers given in figure 39(a). The lOARD input in figure 39(b) 
is identical to the lOARD input presented in figure 38. 
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SUB-STPUCTURES APE DEFINED AS FOLLOWS 
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(b) HCARD and ICARD data used by VIPASA 

Figure 39.- Concluded. 
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